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ABSTRACT 

A  technique  has  been  developed  which  accurately 
predicts  the  performance  of  underwater  acoustic  arrays  prior 
to  array  construction.  The  technique  Is  based  upon  the 
measurement  of  lumped-parameter  equivalent  circuit  values 
for  each  element  In  the  array,  and  Is  accurate  In  predicting 
the  array  transmit,  receive  and  beam  pattern  response. 

The  measurement  procedure  determines  the  shunt 
electrical  and  motional  circuit  elements  from  electrical 
immlttance  measurements.  The  electromechanical 
transformation  ratio  is  derived  from  In-alr  measurements  of 
the  radiating  face  velocity  and  the  Input  current  to  the 
transducer  at  resonance.  The  equivalent  circuit  values  of  a 
group  of  Tonpllz-type  transducers  were  measured,  and  the 
self  and  mutual  Interaction  acoustic  loadings  for  a  specific 
array  geometry  were  calculated.  The  response  of  the 
elements  was  then  predicted  for  water-loaded  array 
conditions. 

Based  on  the  predictions,  a  selection  scheme  was 
developed  which  minimized  the  effects  of  inter-element 
variability  on  array  performance.  The  measured 
transmitting,  receiving  and  beam  pattern  characteristics  of 
a  test  array,  built  using  the  selected  elements,  were 
compared  to  predictions  made  before  the  array  was  built.  The 
results  Indicated  that  the  technique  Is  accurate  over  a  wide 
frequency  range. 
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CHAPTER  1 


INTRODUCTION 


1 . 1  Problem  Definition 

Underwater  acoustic  arrays  are  frequently  comprised  of 
a  number  of  individual  plezoceramlc  transducer  elements, 
each  of  which  is  electrically  Independent.  The  array 
output  Is  the  shaded  (weighted)  sum  of  the  output  of  each 
of  the  elements.  Since  the  transmitting  and  receiving 
characteristics  of  an  array  are  dependent  upon  the 
characteristics  of  the  individual  array  elements,  much  of 
the  design  of  underwater  acoustic  arrays  focusses  on  the 
design  and  testing  of  the  transducer  elements. 

Due  to  a  lack  of  an  adequate  prediction  tool,  the 
design  of  elements  for  an  array  is  often  an  educated  "cut-* 
and-try"  approach.  Involving  the  construction  of  test 
arrays.  The  problem  stems  from  difficulties  encountered  in 
predicting  the  response  of  the  elements  prior  to  array 
assembly.  Current  In-alr  measurement  procedures  do  not 
provide  enough  information  to  generate  model  parameters 
capable  of  the  inclusion  of  a  changeable  acoustic  loading. 
Air  and  water  loading  are  very  different,  and  the  element 
response  varies  considerably  with  the  change  in  medium. 
Further,  the  loading  is  affected  by  mutual  interactions 
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between  elements,  making  element  loading  a  function  of 
position  within  the  array.  Thus  an  element  operating 
singly  will  perform  differently  when  operating  within  an 
array.  Because  of  this,  it  has  not  been  possible  to 
accurately  predict  an  element's  performance  within  an 
array,  prior  to  assembly,  without  construction  of  test 
arrays  . 

Beam  pattern  response  Is  another  area  of  concern  In 
the  array  design  process.  The  array  shading  coefficients 
have  the  effect  of  shaping  the  beam  pattern  of  the  array. 
Varying  the  shading  coefficients  changes  the  beam  pattern. 
The  coefficients  are  altered  by  the  variation  In  response 
of  each  element.  In  order  to  minimize  the  effects  of 
element  selection  and  placement  upon  the  array  shading,  all 
of  the  elements  should  have  similar  response 
characteristics.  However,  due  to  variations  In  assembly 
such  as  glue  joint  thicknesses  and  alignment  problems,  no 
two  elements  will  have  exactly  the  same  response.  These 
response  variations  manifest  themselves  In  the  array  as  a 
degradation  of  the  beam  pattern.  The  use  of  elements  with 
poor  response  uniformity  results  In  Increased  sldelobe 
levels  In  the  beam  pattern. 

Kendlg^  proposed  a  method  for  reducing  the  sldelobe 
levels  In  arrays.  He  devised  a  "scheme”  for  placing  the 
elements  within  the  array  in  order  to  cancel  out  the 
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variations  In  element  performance.  The  "scheme"  Is 
theoretically  sound,  but  the  method  suffers  from  some 
practical  drawbacks,  especially  In  the  prediction  of 
element  performance.  Although  Kendlg  developed  a  means  of 
comparing  elements  to  one  another,  the  comparison  Is  valid 
only  at  the  transducer's  resonances,  and  Is  further 
handicapped  by  a  number  of  simplifications  and  assumptions. 
Even  so,  this  method  has  been  used  to  produce  arrays  with 
lower  sldelobe  le^vels  than  arrays  built  with  a  random 
selection  process. 

A  method  capable  of  predicting  the  response  of 
transducer  elements  within  an  array  is  needed  in  order  to 
overcome  the  drawbacks  in  the  current  design  and  selection 
process;  such  a  method  is  described  In  this  study. 

1 . 2  Goals  of  the  Study 

Using  Interactive  computer  techniques,  this  study 
provides  an  automated  means  for  measuring  transducer 
elements,  predicting  their  responses  when  used  within  an 
array,  and  using  this  information  In  the  Kendlg  selection 
scheme  and  in  the  prediction  of  array  response.  The 
specific  objectives  of  this  thesis  are  : 

1.  Analyze  the  piezoelectric  transducer  element  In  terms 
of  a  lumped  parameter  equivalent  circuit  model, 
capable  of  Including  radiation  loading  effects. 


2.  Develop  a  means  for  automatic  analysis  of  elements  In 
terms  of  the  equivalent  circuit  parameters. 
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3.  Generate  response  characteristics  for  each  transducer, 
having  Included  acoustic  radiation  and  mutual 
Interaction  effects. 

4.  Use  the  response  characteristics  to  find  the  total 
response  for  an  array  of  elements. 

5.  Use  the  response  characteristics  to  Improve  the  Kendlg 
selection  process  for  assigning  array  locations  to  a 
given  set  of  elements. 

6.  Use  the  response  characteristics  to  predict  the  array 
beam  pattern  respoise  for  an  array  of  elements. 


The  first  topic  Is  discussed  In  Chapter  Two;  readers 
already  familiar  with  the  lumped-parameter  equivalent 
circuit  model  may  choose  to  skip  over  this  section. 

Chapter  Three  discusses  the  second  point,  and  the  remaining 
topics  are  covered  In  Chapter  Four. 
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CHAPTER  2 


THE  EQUIVALENT  CIRCUIT  MODEL 


2 . 1  Background 

A  common  type  of  transducer  used  in  underwater 

acoustics  Is  the  "Tonpllz"  -  from  the  German  words  for 

sound  and  mushroom.  Figure  1  shows  the  structure  of  the 

transducer.  Many  equivalent  circuit  model  analyses  have 

2  3 

been  done  using  distributed  parameters,  *  and  are  useful 
In  the  design  of  the  transducer  element.  The  circuit  Is 
valid  over  a  wide  frequency  range,  and  clearly  shows  design 
tradeoffs.  However,  the  circuit  element  parameters  cannot 
be  Inferred  from  tests  of  the  assembled  element,  and  cannot 
account  for  assembly  problems  such  as  misalignment  and 
improper  glue  Joints.  Thus,  the  distributed  parameter 
model  would  not  be  suitable  for  the  purposes  of  this  study. 

The  lumped-parameter  equivalent  circuit  model  Is  valid 
over  a  much  smaller  frequency  range  around  the  transducer's 
mechanical  resonance;  since  It  'lumps"  the  effects  of  the 
Individual  transducer  components  together.  It  Is  somewhat 
less  useful  In  the  design  phase.  However,  the  equivalent 
circuit  element  parameters  can  be  derived  from  tests  of  the 
assembled  transducer.  Further,  the  transducers  are 
commonly  used  near  or  below  their  fundamental  resonance 
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frequency,  where  the  equivalent  circuit  Is  valid.  For 
those  reasons,  the  lumped  parameter  equivalent  circuit 
model  will  be  used. 

A  common  starting  point  In  the  equivalent  circuit 
analysis^  Is  a  simple  block  of  piezoelectric  material,  one 
quarter  wavelength  thick,  one  face  rigidly  mounted  and  the 
other  radiating  Into  the  medium.  The  analysis  leads  to  an 
equivalent  circuit,  but  the  values  of  the  parameters 
derived  are  not  suitable  to  the  Tonpllz  configuration. 
Cherpak^  has  done  an  analysis  of  different  Tonpllz  type 
transducers,  but  again,  the  results  neither  predict  nor 
account  for  the  element  variability  found  In  actual 
practice. 

Our  purpose,  then.  Is  to  derive  the  lumped-parameter 
equivalent  circuit,  and  to  put  the  values  of  the  parameters 
In  terms  of  measurements  conducted  on  assembled  units. 

From  the  references  cited,  we  can  be  confident  that  the 
circuit  Is  topologically  correct;  tests  will  later  show  the 
validity  of  the  model  In  simulating  actual  element 
per  f o  rmanc e  . 

2 . 2  Assumptions  and  Limitations 

In  order  to  derive  a  model  of  the  transducer.  It  Is 
necessary  to  make  some  assumptions  about  the  transducer's 
behavior  and  structure.  Attempting  to  model  the  transducer 
completely  would  lead  to  an  unnecessarily  complex  and 
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unwieldy  equivalent  circuit.  Further,  the  model  circuit 
elements  are  to  be  derived  from  measurements  of  actual 
transducers.  The  number  of  suitable  measurements  poses  a 
limit  on  the  number  of  circuit  elements  In  a  model  that  can 
be  adequately  derived.  Therefore,  certain  assumptions  and 
restrictions  will  be  made  on  the  model.  Imposing 
restrictions  generates  limitations  on  the  model's  accuracy 
and  suitability  to  the  array  design  process. 

The  following  restrictions  apply  to  the  model  of  the 
element's  performance: 

1.  The  entire  circuit  will  be  linear,  passive,  and 
rec Iprocal  . 

2.  The  mechanical  domain  can  be  treated  as  a  simple 
mechanical  oscillator  (one  degree-of-freedom)  . 

3.  The  radiating  face  of  the  transducer  has  a  fixed 
velocity  distribution  (one  degree-of-freedom). 

4.  The  electrical  domain  has  little  dissipation. 

The  first  restriction  limits  the  possible  model 
elements  to  pure  resistances,  capacitances,  and  Inductances 
In  the  electrical  domain,  and  to  pure  springs,  masses, 
etc..  In  the  mechanical  and  acoustical  domains.  Although 
It  will  become  apparent  that  linear  elements  are  not  able 
to  completely  model  the  transducer's  behavior  over  a  wide 
frequency  range,  they  are  sufficient  over  the  range  of 
Interest,  up  to  the  fundamental  resonance.  By  limiting  the 
analysis  to  linear  elements,  the  transducer's  behavior  can 
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be  described  by  sets  of  linear  differential  equations.  For 
the  harmonic  steady  state,  the  linear  differential 
equations  may  be  transformed  Into  linear  algebraic 
equations  In  the  frequency  domain.  Because  of  this,  the 
circuit  variables  will  be  complex  phasors,  and  the  the 
circuit  element  values  will  be  complex  Impedances. 

The  second  restriction  limits  the  model  to  a  single 
mechanical  resonance.*  However,  since  transducers  are 
typically  operated  only  below  or  near  that  resonance, 
providing  for  more  degrees-of-freedom  and  resonances  does 
not  significantly  Increase  the  model's  accuracy.  The 
additional  degrees-of-freedom  would  account  for  such  things 
as  higher  harmonics  of  the  fundamental,  transverse 
vibrations,  and  spurious  resonances.  These  generally  are 
not  of  great  concern  nor  are  they  easily  measured. 

The  third  assumption  requires  that  the  transducer's 
radiating  face  moves  as  a  rigid  piston.  Therefore  this 
model  will  not  treat  a  transducer  with  head  flap.^ 

The  final  restriction  Is  true  for  transducers  of 
Interest,  constructed  of  plezoceramlc  material.  One 
problem  with  Including  electrical  losses  Is  that  they  are 
hysteretlc,  and  therefore  nonlinear.  If  linearized  over 
some  frequency  range,  the  loss  turns  out  to  be  small  enough 

*It  will  be  shown  how  additional  degrees-of-freedom  and 
resonances  can  be  added. 
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to  neglect  In  comparison  to  other  circuit  elements.  The 
Implications  of  this  restriction  will  be  discussed  further 
In  the  sections  dealing  with  measurement  procedures. 

A  generalized  linear  equivalent  circuit  model  of  the 
transducer  can  be  generated  from  a  simple  consideration  of 
the  transducer's  function.  The  function  of  a  transducer  Is 
the  conversion  of  energy  from  the  acoustical  domain  to  the 
electrical,  and  vice  versa.  In  the  case  of  piezoelectric 
transducers,  acoustical  energy  Is  first  converted  to 
mechanical  energy;  mechanical  energy  Is  then  converted  Into 
electrical  energy.  It  Is  necessary  for  the  lumped- 
parameter  equivalent  circuit  model  of  a  piezoelectric 
transducer  to  account  for  the  transducer's  properties  In 
the  electrical,  mechanical,  and  acoustic  domains  ,  as  well 
as  the  transformations  between  these  domains.  With  this  In 
mind.  Figure  2  Is  a  generalized  model  of  a  transducer 
element.  There  are  five  two-port  networks  to  represent  the 
electrical,  mechanical,  and  acoustical  domains,  and  the  two 
transformations  between  the  domains.  Note  the  different 
potential  and  flux  variables  used  between  the  networks. 

This  model  will  be  made  more  specific  In  the  next  section, 
using  the  above  restrictions. 

2 . 3  Derivation  of  Model  Circuit  Elements 

Starting  with  the  generalized  linear  equivalent 
circuit  model  In  Figure  2,  the  exact  structure  of  each  of 
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the  networks  can  be  determined  from  the  above  restrictions 
and  an  Intuitive  look  at  the  transducer's  construction. 

The  structure  of  the  first  network,  representing  the 
electrical  domain,  can  be  deduced  from  the  physical 
structure  of  the  transducer,  Figure  1.  Neglecting  the 
piezoelectric  properties  of  the  ceramic,  the  electrodes 
form  a  simple  parallel  plate  capacitor  stack,  with  the 
ceramic  serving  as  a  dielectric  material.  The  capacitance 
Is  a  function  of  the  number  and  area  of  the  electrodes,  and 
the  thickness  and  dielectric  constant  of  the  ceramic 
sections.  There  Is  also  a  small  dielectric  loss  from  the 
conductance  of  the  ceramic.  The  network  representing  the 
electrical  domain  consists  of  a  capacitor  In  parallel  with 
a  resistor. 

The  capacitive  property  arises  from  the  ability  of  a 
dielectric  material  to  hold  a  potential  charge  across  a 
pair  of  electrodes.  The  ceramic  In  a  transducer,  however, 
is  both  dielectric  and  piezoelectric.  The  charge  across 
the  ceramic  therefore  also  causes  a  mechanical  deformation 
In  the  material's  crystalline  structure,  coupling  the 
electrical  and  mehanlcal  domains. 

The  equivalent  circuit  model  for  the  electromechanical 
transformation  Is  found  by  examining  the  nature  of  the 
piezoelectric  phenomenon.  Mechanically  deforming  a 
piezoelectric  material  causes  a  charge  separation  within 
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Its  crystalline  structure.  The  voltage  potential  caused  by 
the  separation  will  vary  with  the  orientation  of  the 
mechanical  stress  and  voltage  measurement  with  respect  to 
the  crystal  polarization  orientation.  The  voltage  Is  a 
function  of  the  size  of  the  material,  the  amount  of 
deformation,  and  the  piezoelectric  stress  factor 
appropriate  to  the  crystal  polarization  orientation.  By 
Including  the  compliance  coefficient  of  the  plezoceramlc 
material  Into  the  expression,  the  transformation  factor 
takes  on  the  dimension  of  f orce/voltage ,  In  units  of 
Newtons/Volt.  By  further  substitution,  the  ratio  Is 
equivalent  to  current/ velocity  (Amp/Meter/Second) .  The 
equivalent  circuit  element  that  expresses  this  type  of 
transformation  Is  an  Ideal  transformer.  Figure  3  shows  the 
Ideal  transformer  and  the  circuit  relationships  applicable 
to  It.  By  convention,  the  electromechanical  transformation 
factor  Is  given  the  symbol  ip  ,  (A/M/Sec  or  N/V),  and  the 
turns  ratio  Is  expressed  as  1/<I>:1  or  1  :  4* ,  which  Is 
equivalent  to  letting  C  In  Figure  3  equal  l/4>. 

From  the  second  restriction,  the  network  representing 
the  mechanical  domain  Is  a  simple  mechanical  oscillator. 

The  representation  of  a  mechanical  oscillator  Is  the  system 
shown  In  Figure  4A,  consisting  of  a  mass,  a  spring,  and  a 
dissipative  element.  Under  the  assumption  that  the 
dissipation  Is  small,  and  Is  proportional  to  the  velocity, 
the  equivalent  circuit  for  the  mechanical  system  Is  shown 


B  C:1  D 


TRANSFORMATION  RELATIONSHIPS 
POTENTIAL:  G  =  A/C 
FLUX:  D  =  C/B 

IMPEDANCE:  G/D  = 

STEP  DOWN:  C>1  ^ 

STEP  UP:  C  <1 


Figure  3.  The  Ideal  Transformer 


Figure  4A.  A  Mechanical  Oscillator 

Figure  4B.  The  Electrical  Equivalent  of  a  Mechanical 

Oscillator 
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In  Figure  4B.  The  circuit  elements  are  In  series  because 
the  velocity  of  each  element,  measured  at  point  A,  Is  the 
same.  Since  velocity  is  the  flux  variable  In  the  model, 
the  same  velocity  must  “flow  through”  each  of  the 
equivalent  circuit  elements.  By  connecting  In  parallel 
sections  like  the  one  shown,  additional  degrees-of-f reedom 
and  resonances  may  be  added. 

Anticipating  the  results  of  the  next  paragraph,  the 
mechanoacoustlcal  transformation  Is  accomplished  by  the 
radiating  face  of  the  transducer.  Under  the  assumption  of 
"lumped"  elements,  and  looking  at  Figure  1,  the 
displacement  generated  by  the  ceramic  elements  must  be 
equal  to  the  displacement  of  the  radiating  face,  taken  with 
respect  to  the  tail  mass.  must  equal  U^^ .  Therefore, 

the  equivalent  circuit  for  the  mechanical  domain  In  the 
model  Is  shown  In  Figure  5. 


The  mechanoacoustic  transformation  Is  accomplished 
through  the  area  of  the  radiating  face  of  the  element.  A 
displacement  of  the  face  causes  a  change  In  volume  In  the 
acoustic  medium:  the  velocity  of  the  face  produces  a 

volume  velocity  In  the  medium.  Similarly,  a  force  F^^g 
on  the  face  becomes  a  force  per  unit  area,  or  pressure  Pg 
In  the  medium.  Mathematically,  this  Is  expressed  as 


SU 


ma 


and 


(1) 
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where  S  Is  the  area  of  the  face.  The  equivalent  circuit 

representation  for  the  mechanoacoustlc  transformation  Is  an 

Ideal  transformer  with  a  turns  ratio  of  S:l,  and  S  has  the 
2 

units  of  Meters  . 

The  final  network,  representing  the  acoustical  domain, 
cannot  be  strictly  analysed  as  a  combination  of  lumped 
circuit  elements.  The  solution  of  the  radiation  loading  on 
a  flat  piston  Is  very  complex,  and  Is  further  complicated 
by  mutual  coupling  effects  In  arrays.  However,  for  a  given 
frequency  and  array  velocity  distribution  (shading 
pattern),  there  Is  a  unique  radiation  loading  on  each 
element*,  which  may  be  represented  as  a  generalized 
Impedance  .  Figure  6  Is  the  circuit  representation  for 
the  acoustical  domain.  The  power  output  of  the  transducer 
Is  proportional  to  the  volume  velocity  through  the  real 
part  of  Zj^,  labeled  . 

Figure  7  then.  Is  the  lumped-parameter  equivalent 
circuit  model  of  a  Tonpllz  piezoelectric  transducer. 

2 . 4  Input/Output  Equations  for  the  Model 

The  lumped-parameter  equivalent  circuit  model  can  be 
analysed  by  standard  electrical  network  techniques  to  find 
the  Input/output  equations  of  the  transducer.  Figure  8 


*The  techniques  for  finding  the  radiation  loading  and 
mutual  coupling  Impedance  are  discussed  In  Appendices  A  and 
B,  respectively. 


nsmlt  and  Receive  Cases 
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results  from  removing  the  transformers  by  the  appropriate 
multiplication.  Figure  8A  Is  the  circuit  for  the  transmit 
case  and  Figure  8B  is  for  the  receive  case. 

For  the  receive  case,  the  measured  electrical  quantity 
can  either  be  the  open  circuit  voltage  8^^^,  or  the  short 
circuit  current  Similarly,  In  the  transmit  case,  the 

driving  electrical  source  can  be  modelled  either  as  an 
Ideal  voltage  source  or  an  Ideal  current  source 

Although  the  acoustical  domain  cannot  be  considered  either 
shorted  or  opened,  the  transmit  and  receive  cases  use 
different  circuit  parameters.  In  the  receive  case,  the 
Input  variable  la  the  pressure  Pin  at  the  face  of  the 
transducer;  In  the  transmit  case,  the  output  variable  Is 
the  volume  velocity  produced  by  the  transducer. 

By  convention,  for  the  receive  case,  the  output 
voltage  and  current  sensitivities  for  the  transducer 

element  may  be  expressed  as 

'Pe'k  -  lE^ut'/IPlnl  'Ri'k-  Hout'/lPlnl  <2) 

Similarly,  the  transmit  response  may  be  written 

I  Te  I  k  ■  I '^out  ^  ^  ^  ®in  ^  ^^l^k  “  ^ '^out  ^  ^  ^  ^  in  ^ 

Since  the  derived  equivalent  circuit  Is  reciprocal  and 
passive.  It  can  be  shown  that^ 

^*^i^k  "  ^^e^k  ^  ^e  ^  k  “  ^^l^k 


(4) 


Note  that  this  Is  not  the  same  as  acoustic  reciprocity, 
which  Includes  a  factor  accounting  for  sound  transmission 
Into  the  medium. 


Solving  the  equivalent  circuit  for  the  first  pair 
yields 


IR 


tiiC  6S 

I  _ _ 1_ _ 

i  R  9  99  99  99% 

[(l-U)^C  M  -  wC  S  X  )  +  t/c  (R  +  S  R  )  ] 

mm  ma  m  m  a' 


(5) 


,  ojC  (R  +  S  R  ) 

Arg(Ri)^  “  tan"^  { - l-Ji! - ® - —  } 

1  -  oj^C  M  -  U)C  S'^X 
mm  m  a 

If  the  transducer  is  operated  at  a  frequency  such  that 


-  1  equals  zero,  then  the  equations  reduce  to 


iRlIk  ■  <0S/[r2  +  ZR^R^S^  + 

Arg(Ri)k  •  tan‘M-[Rnj  +  S^R^J/S^x^} 


(6) 


where  I I  »  [R^  +  X^]^.  Solving  the  circuit  for  (Re)^ 
(T^)k  yields 


IR 


(j)S 


e'k’oo  77  7  22^ 

[((t*  -  M  +  C  /C  -  cue  S  X  )  +  toe  (R  +S  R,)  ] 

‘  om  om  oa  om  a  ^ 

,  U)C  (R  +  S^R  ) 

Arg(Rg)k  •  tan-M^ - 2_E - 2 -  } 

(t)  -  0)  C  M  +  C  /C  -  wC  S  X 

^  o  m  o  m  o  a 

If  the  transducer  Is  operated  at  the  frequency  such 


that  1  -  +  Cg/Cju'^^  equals  zero,  then  the 


-  0)^ 

equations  simplify  to 


iRglk  ■  <i)S/(u)C„)[R2  +  2R„R^s2  + 
Arg(Re)k  '  tan'M-fR^  +  RgS^j/s^x^} 


(8) 
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These  are  the  input/output  or  transfer  function 
equations  for  the  transducer.  Equations  5  and  7  will  be 
used  for  predicting  the  amplitude  and  phase  of  the 
transducers  in  the  array  design  process.  The  choice  of  the 
frequencies  which  reduced  Equations  5  and  7  are  not 
arbitrary:  the  frequencies  will  be  shown  to  be  the 
resonant  and  antlresonant  frequencies,  respectively. 

2 . 5  The  Reduced  Model 

For  measurements  of  the  electrical  Immlttance 
properties  of  the  transducer,  the  complete  equivalent 
circuit  model  may  be  reduced  to  a  more  convenient  form.  If 
we  assume  that  the  transducer  is  measured  in  air,  for  which 
the  radiation  loads  are  relatively  small,  then  R  and  X 
can  be  set  to  zero.  Similarly,  we  can  assume  that  there  Is 
no  significant  pressure  field  at  the  transducer  face. 
Applying  those  assumptions  to  Figure  SB,  we  can  rename  the 
parameters  as  traditional  electrical  circuit  elements.  The 
results  are  shown  In  Figure  9,  where 

^o  "  ^o  ^  ^ 

This  simple  circuit  Is  sufficient  to  describe  the 
electrical  Immlttance  behavior  of  the  transducer,  with  the 
assumption  that  and  can  be  neglected.* 

*  This  was  tested  as  part  of  the  experiments  described  In 
Section  3.2. 


Figure  9.  The  Reduced  Equivalent  Circuit  Model 
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It  will  be  showa  that  all  of  the  non-acoustic 
parameters  In  the  complete  equivalent  circuit,  except  the 
transformation  factor  (|)  ,  can  be  found  by  treating  the 
transducer  as  a  one-port  device.  The  technique  for  finding 
<p  Is  discussed  In  Section  3.3. 

2 . 6  One-Port  Behavior  of  the  Reduced  Model 

The  reduced  equivalent  circuit  model  can  be  used  to 
predict  the  electrical  Immlttance  response  of  the 
transducer  over  a  frequency  band  around  the  resonance 
frequency.  The  values  of  the  circuit  elements  can  be  found 
from  measurements  of  the  transducer's  electrical 
properties,  and  relating  the  measurements  to  the 
predictions  of  the  model.  The  following  discussion  will  be 
limited  to  the  reduced  equivalent  circuit  model,  as  shown 
In  Figure  9. 


From  the  figure,  the  Input  electrical  admittance  of 
the  transducer  Is: 


Y  m  itiiC  +  J  (i)C  (1  —  LCcu  )  ^  fjQ  C  R 
J  O  ^  2  ?  2  2  2 

(i  -  LCu  )  +  U)  C  R 


(10) 


This  will  show  two  resonances  :  one  from  the  motional 
branch  alone,  and  one  from  the  coupled  system.  The 
motional  resonance  frequency  Is 

% 


OJy  -  l/(LC) 


(11) 


At  this  frequency,  the  admittance  Is 
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*Ly  ■  JVo  +  <“> 

The  second,  or  antiresonance  frequency  Is 

-  [(C  +  Co)/(LCCo)]^  (13)- 

The  conductance  of  the  transdcuer  Is  a  maximum  at  the 

resistance  of  the  transducer  Is  approximately  a  maximum  at 

Note  that  Eqs.  11  and  13  are  the  same  as  the 
conditions  used  to  reduce  Eqs.  5  and  7. 

Figure  10  shows  the  Input  admittance  versus  frequency 
for  Che  model  near  Che  motional  resonance.  The  susceptance 
curve  Is  Che  sum  of  the  admittance  of  the  capacitor  (the 
"blocked”  susceptance)  and  the  admittance  of  Che  motional 
branch.  The  conductance  curve  Is  due  to  the  motional 
branch  alone.  The  peak  In  the  conductance  curve  occurs  at 
the  resonance  frequency,  At  this  frequency,  the 

motional  susceptance  goes  to  zero;  from  Eq.  12,  the 
motional  conductance  Is  just  Che  Inverse  of  Che  motional 
resistance  element  R. 

The  two  frequencies  labeled  and  ^2  half¬ 
power  or  quadrantal  frequencies.  Together  with  the 

points  define  the  or  quality  factor  of  the  motional 

branch,  which  Is  expressed  as 

-  u)y/(iui  -  u)2)  - 
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From  basic  circuit  theory,  the  may  also  be  expressed  In 
terms  of  the  motional  branch  circuit  elements: 

■  l/UyRC  (15) 

In  the  absence  of  loading,  the  Is  a  measure  of  the 
mechanical  dissipation  of  the  transducer  at  resonance.  A 
higher  Implies  a  more  efficient  motional  branch. 

Another  measure  of  a  transducer's  performance  Is  the 
ratio  of  the  energy  stored  In  the  motional  branch  to  the 
total  energy  stored  In  the  transducer.  By  examining  the 
model  at  low  frequencies,  where  the  capacitive  elements 
dominate  the  Impedance,  the  electromechanical  coupling 
coefficient,  called  Rgff*  becomes  obvious.  From  Figure  11, 
the  energy  In  the  motional  branch  Is  (1/2)CV  ,  and  the 
total  energy  stored  Is  (1/2)(C  +  .  Therefore, 

kgff  -  [C/(C  +  -  [1  -  (16) 

The  above  analysis  has  shown  that  the  reduced 
equivalent  circuit  model  has  two  resonances,  at  Uy  and 
From  those  frequencies  and  the  half-power  points,  the 
and  kg££  are  derived.  Section  3.1  will  discuss  the 
experimental  techniques  for  measuring  transducers  to  find 


these  parameters 


ELECTRICAL  MOTIONAL 

BRANCH  BRANCH 


Figure  11.  Low  Frequency  Reduced  Mode 


CHAPTER  3 


THE  MEASUREMENT  OF  CIRCUIT  PARAMETERS 


3 . 1  Input  Immlttance  Measurements 

In  the  previous  section,  the  Input  Immlttance  behavior 
of  the  reduced  equivalent  circuit  model  was  derived.  Using 
Equation  10  as  a  starting  point,  the  resonance  and 
antiresonance  frequencies,  the  and  the  coupling 

coefficient  were  derived  In  Equations  11  through  16.  The 
purpose  of  this  section  Is  the  derivation  of  the  equations 
needed  to  find  the  equivalent  circuit  values  from 
measurements  of  Input  Immlttance. 

From  Equation  12,  If  By  and  Gy  are  the  susceptance  and 
conductance  of  the  circuit  at  Wy,  then 

R  •  1/Gy  and  C^  -  By/cOy  (17) 

From  the  measurement  of  G  and  Q  , 

y  ^m 

Finally,  the  value  of  the  Inductor,  L,  is  derived  from  C 
and  <*)y, 


L  -  1/U)^C 


(19) 
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These  four  equations  are  sufficient  for  analysis  of 
the  equivalent  circuit  model. 

In  practice,  the  circuit  element  values  are  found  by  a 
slightly  different  set  of  relations,  for  reasons  of 
repeatability  and  ease  of  measurement.  If  the  transducer 
element  Is  measured  at  a  low  enough  frequency,  then  the 
Impedance  Is  dominated  by  the  capacitances.  The  low 
frequency  susceptance  Is  due  to  the  combined  effect  of  C 
and  Cq .  Therefore, 

C  +  Cq  »  B/(i)  aj<<  u)y  (20) 

This  equation  Is  used  with  Eq .  18  to  find  ,  the  shunt,  or 
clamped  capacitance. 

Figure  12  Is  the  equipment  set-up  used  for  the 
Impedance  measurement.  The  HP3750  Impedance  Analyzer  Is 
designed  to  make  Insertion  loss/gain  and  phase  angle 
measurements  on  two-port  devices.  It  was  used  to  find  the 
Input  Immlttance  of  one-port  devices  by  use  of  a  special 
Input  termination  and  use  of  the  appropriate  formulae 
stored  In  the  calculator.  The  details  of  the  termination 
and  formulae  are  discussed  In  Appendix  C.  The  HP9825B 
calculator  was  programmed  to  measure  the  admittance  and 
Impedance  characteristics  of  the  transducer  element  under 
test,  and  to  calculate  the  equivalent  circuit  values  for 
the  element.  The  low  frequency  capacitance  measurements 
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were  taken  at  three  frequencies  near  00^/20,  and  then 
averaged.  The  frequencies  u  and  o)  were  found  by  sweeping 
over  a  specified  frequency  range,  until  the  frequency  of 
maximum  conductance  and  the  frequency  of  maximum 

resistance  (<^2^  were  located.  The  half~power  frequencies 
were  specified  as  the  frequencies  for  which  the  conductance 
was  one-half  G^.  From  these  measurements,  the  calculator 
used  Equations  17  through  20  to  find  the  R,  L,  C,  and  of 
the  transducer  element.  The  circuit  element  values  were 
stored  In  the  calculator  memory  for  later  use.  Plots 
showing  the  magnitude  and  phase  of  the  Impedance  versus 
frequency  could  also  be  made.  The  solid  curve  In  Figure  13 
Is  an  example  of  a  typical  measured  Impedance  plot. 

To  establish  the  reliability  and  repeatability  of  the 
measurements,  many  tests  were  repeated  on  a  control  group 
of  transducer  elements.  Once  mounting  and  measurement 
techniques  were  refined,  a  group  of  over  60  transducer 
elements  were  tested  and  the  results  stored  for  later  use. 

3 . 2  Verification  of  the  Immlttance  Techniques 

In  order  to  verify  the  Immlttance  measurements,  a 
number  of  calibrations  and  built  In  self-checks  are 
calculated . 


The  accuracy  of  the  low  frequency  capacitance  was 
verified  using  the  Impedance  program  to  measure  a  number  of 


and  Measured  Transducer  Impedance 


il, 
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capacitors,  whose  values  were  In  the  range  of  the 
transducer  elements  tested.  The  readings  were  found  to  be 
within  +.6  percent  of  capacitance  bridge  measurements. 

It  was  assumed  In  Section  2.5  that  the  acoustic 
radiation  loading  In  air  Is  sufficiently  low  that  R_  and  X. 

a  cl 

could  be  set  to  zero.  This  was  tested  by  placing  elements 
In  a  bell  jar,  and  measuring  the  Immlttance  of  the  elements 
under  atmospheric  and  low  pressure  conditions.  The  results 
Indicated  that  air  loading  caused  a  .5  percent  error  In  the 
measurement  of  R.  There  was  no  significant  shift  In  the 
measured  values  of  any  of  the  other  parameters  of  Interest. 
It  was  therefore  concluded  that  the  assumption  of 
negligible  In-alr  loading  was  correct. 

A  self-check  routine  In  the  program  calculated  and 

printed  the  Kgff  of  the  transducer  using  both  terms  of 

Eq.  16.  By  comparing  the  value  from  the  capacitance 

measurements  and  the  value  from  the  u)  ,  and  uo  measurements, 

y  z  ’ 

the  self-consistency  of  the  entire  measurement  was  checked. 
The  differences  were  typically  less  than  .7  percent  for 
most  of  the  elements.  This  Indicated  the  resonance 
frequency  and  the  capacitance  measurements  were  Indeed 
consistent . 

Some  of  the  elements  had  a  portion  of  the  head  mass 
removed  to  allow  for  better  array  geometry.  These  elements 
also  had  their  tall  masses  changed  or  "trimmed"  to 


compensate  for  Che  shift  in  ui^  caused  by  Che  loss  of  head 
mass.  These  changes  In  mass  caused  Che  measured  and 
values  Co  shift,  and  therefore  changed  the  calculated  value 
of  C  for  the  transducers.  When  the  values  of 
calculated  from  C  and  were  compared  to  values  calculated 
from  (jLi„  and  (ii  ,  there  was  a  discrepancy  of  7  percent.  An 
examination  of  Che  measured  impedance  plots  revealed  that 
these  elements  have  an  extra  resonance  slightly  below  , 
thus  affecting  their  behavior  near  resonance.  The 
equivalent  circuit  model  as  derived  cannot  account  for  the 
behavior  of  these  elements,  because  of  the  extra  resonance. 
Therefore  caution  was  taken  when  interpreting  the  model's 
response  prediction  for  these  elements. 

A  second  self-check  routine  had  two  complementary 

purposes:  to  test  the  accuracy  of  the  measurement 

procedure,  and  to  test  the  accuracy  of  the  equivalent 

circuit  model.  To  do  this,  the  calculator  was  programmed 

using  Equation  10  with  the  measured  equivalent  circuit 

values  Inserted.  The  immlttance  curves  were  redrawn  using 

the  lumped  paramenter  equivalent  circuit.  The  results  were 

both  printed  and  plotted.  The  printed  results  showed  that 

the  technique  could  find  the  equivalent  circuit  element 

values  to  within  .5  percent.  The  plots  showed  the  the 

equivalent  circuit  accurately  modeled  the  transducer's 

input  Impedance  over  a  range  near  ,  as  seen  in  Figure  13. 

The  model  was  not  as  accurate  in  the  range  near  cj  , 

z 


where 
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the  predicted  Impedance  was  generally  30  percent  higher 
than  measured.  This  discrepancy  arises  from  slight 
Inaccuracies  In  the  measurement,  and  In  the  Inaccuracies 
caused  by  neglecting  R^.  Appendix  D  discusses  the  error 
caused  by  the  removal  of  R^,  and  suggests  a  technique  for 
reducing  It.  It  was  found,  however,  that  this  error  had 
little  effect  on  the  final  results  of  the  study. 

3 . 3  Velocity-Current  Measurements 

The  electromechanical  transformation  factor  cannot  be 
found  by  measurement  of  Input  electrical  properties  alone. 
Since  the  Immlttance  measurement  treats  the  transducer  as  a 
one-port  device  (Fig.  9),  It  cannot  get  at  "both  sides"  of 
the  electromechanical  transformation.  Therefore,  It  Is 
necessary  to  measure  some  electrical  and  mechanical 
properties  of  the  transducer  simultaneously. 


Referring  to  Figure  7,  the  electromechanical 
transformation  factor  can  be  expressed  either  as  a 
voltage/f  orce  (Eg^j/F^j^)  or  a  current/velocl  ty 
ratio.  Although  the  Input  voltage  can  be  easily  measured, 
the  output  force  from  the  Ideal  transformer  cannot; 
however,  the  velocity  of  the  radiating  face  and  the  Input 
current  can  be  measured  simultaneously.  The  current  Ig|||  In 
the  standard  circuit  model  Is  the  same  as  the  motional 
current  In  the  the  reduced  circuit  model.  The  reduced 
circuit  model  can  be  solved  for  the  motional  current 

m 


using  the  element  values  found  by  Impedance  measurements* 
The  motional  current.  In  terms  of  the  electrical  Input 
current  I ^ ,  Is  : 

ll'l  -  ll^l  /  t(l  -  +  C^,/C)2  +  (coCoR)^]^ 

(21) 

Arg(l')  -  Arg(I^)  +  tan~^  f -wC^R/ { 1  -  +  C^/C)] 

The  current  U„  out  of  the  electromechanical 
m 

transformer  equals  the  current  flowing  Into  the 

mechanoacous t Ic  transformer.  This  last  transformer 
represents  the  transducer's  radiating  face  ;  on  one  side  Is 
the  velocity  of  the  mechanical  domain,  and  on  the  other  Is 
the  volume  velocity  of  the  acoustical  domain.  The  two  are 
related  by  the  area  of  the  face.  In  this  case,  only  the 
velocity  of  the  face,  U^,  Is  of  Interest  and  not  the  volume 
velocity  produced  by  the  face.  The  current  can  be 

found  by  measuring  the  velocity  of  the  face.  If  Is  the 
measured  face  velocity,  then  the  electromechanical 
transformation  factor  Is 

<|)  ■  l'  /  Uf  (22) 

The  face  velocity  Is  determined  from  measurements  of 
the  face  displacement  amplitude.  In  the  steady  state,  the 
two  are  related  by  the  radian  frequency.  A  non-contacting 
fiber  optic  probe  Is  used  to  measure  the  displacement 

Q 

amplitude.  Figure  14  Illustrates  the  principle  of 


41 

operation  of  the  light  fiber  probe.  As  a  test  surface 
moves  away  from  the  probe  tip,  the  area  of  overlap  between 
the  transmitted  light  cone  and  the  received  light  cone 
Increases.  The  sensor  electronics  convert  the  light 
amplitude  signal  from  the  probe  Into  a  voltage  proportional 
to  displacement.  By  differentiating  the  displacement 
signal,  the  velocity  of  the  surface  Is  measured. 

Figure  15  shows  the  experimental  set-up  used  to 
measure  the  transformation  ratio.  The  transducers  were 
mounted  by  their  stress  bolts;  Immlttance  measurements  had 
shown  that  the  method  did  not  significantly  alter  the 
transducer's  resonances.  This  mounting  method  also  Insured 
that  the  face  displacement  was  the  sum  of  the  expansions  of 
the  ceramic  sections,  assuming  that  little  of  the 
displacement  Is  taken  up  in  compliant  deformation  of  the 
head  mass,  the  tall  mass,  and  the  stress  bolt. 

The  Input  current  to  the  transducer  was  measured  by  a 
current-to-vol tage  conversion  coll.  The  unit  has  a 
sensitivity  of  .1  volt/ampere,  and  has  sufficient  frequency 
response  to  cover  the  range  of  Interest. 

The  fiber  optic  probe  was  calibrated  by  using  a 
calibrated  positioner  and  a  voltmeter.  The  probe  tip  was 
manually  displaced  In  steps  away  from  the  transducer  face, 
and  simultaneously  the  voltage  from  the  sensor  unit  was 
recorded  at  each  step.  This  provided  a 
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vol tage/dlsplacement  curve  which  was  matched  to  the 
reflectivity  of  the  transducer  face.  It  also  clearly 
showed  the  probe's  range  of  linear  operation. 

During  measurements,  the  output  voltage  from  the 
sensor  unit  Is  amplified  and  filtered.  The  amplification 
raises  the  signal  to  well  above  the  Instrumentation  noise 
floor;  the  filtering  removes  both  the  DC  bias  arising  from 
the  static  offset  of  the  probe,  and  the  high  frequency 
leakage  from  the  sensor  unit  electronics. 

If  Ej  Is  the  displacement  voltage  measured  In 
dB  re  .IV,  A  Is  the  amplification  In  dB,  Is  the  filter 
loss,  S  Is  the  sensitivity  of  the  probe  In  meters/volt,  and 
u)  Is  the  radian  frequency,  then  the  velocity  of  the  face  Is 

Uf  -  (jJS[10(^d  -  A  +  Lf  -  20)/20j  (23) 

This  equation  was  programmed  Into  the  calculator  to  find 
the  face  velocity. 

The  results  of  the  tests  showed  that  the  measurement 
procedure  Is  repeatable  to  within  +2  percent.  The  error 
comes  from  the  repeatability  of  the  manual  positioner 
calibration.  The  error  could  be  reduced  by  designing  a 
more  accurate  and  repeatable  positioning  system.  The 
velocity  proportional  and  current  signals,  as  displayed  on 
an  oscilloscope,  were  free  from  noise  or  distortion.  The 
average  value  of  the  transformation  factor  was  .649  Amp- 
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Sec/Meter;  the  maximum  was  .822  and  the  minimum  was  .537  A- 
S/M. 


3 . 4  Frequency  Dependence  of  the  Transformation  Factor 

In  order  to  verify  some  of  the  assumptions  about  the 
velocity-current  technique,  a  test  was  conducted  to  measure 
the  frequency  dependence  of  the  measured  <P  values.  The 
calculator  was  programmed  to  store  the  measured  velocity 
and  current  of  a  transducer  at  200  frequency  points  over 
a  given  range  around  the  transducer's  resonance.  The 
results  were  corrected  by  using  Equation  21,  and  are  shown 
In  Figure  16.  Beyond  the  frequency  range  shown,  the 
velocity  signal  becomes  significantly  degraded  by  system 
noise. 


The  overall  non-zero  slope  of  the  <P  curve  Is  probably 
due  to  small  Inaccuracies  In  the  measurement  of  R  In  the 
Immlttance  technique.  Since  the  velocity-current 
measurement  Is  done  near  resonance,  the  term  of  Eq .  21 

dominates  the  correction  expression.  Small  changes  In  R 
affect  the  correction  term  and  cause  a  change  In  the  slope 
of  the  <j)  curve.  The  discrepancies  near  Fy  are  probably  due 
to  problems  In  the  measurement  equipment  which  could  not 
keep  pace  with  abrupt  changes  In  current  and  velocity.  The 
of  the  transducer  In  air  Is  high,  and,  near  resonance, 
the  Internal  transducer  parameters  are  changing  too  rapidly 
for  the  measurement  technique  programmed  Into  the 


calculator 


u 


FREQUENCY,  F/Fy  (kHz) 


Figure  16.  Electromechanical  Transformation  Factor  as 

Function  of  Frequency 
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Aside  from  those  discrepancies  Just  mentioned,  the 
curve  Indicates  that  the  electromechanical  transformation 
factor  Is  constant  with  frequency,  as  assumed  In  the 
derivation  of  the  equivalent  circuit  model. 


CHAPTER  4 


APPLICATIONS  TO  ARRAY  PERFORMANCE 


4 . 1  The  Inclusion  of  Radiation  Loading 

The  characteristic  responses  of  a  transducer  element 
In  water  are  vastly  different  than  In  air:  e.g.,  In  water, 
the  resonance  frequencies  shift,  and  the  quality  factor, 

Qjjj,  goes  down.  The  changes  are  caused  by  the  dramatic 
difference  In  the  radiation  loading,  or  acoustic  Impedance, 
between  the  two  media. 


The  radiation  loading  may  be  expressed  as  follows: 


where  p  Is  the  density  of  the  medium,  c  Is  the  speed  of 
sound  In  the  medium,  A  Is  the  area  of  the  radiator,  Rgg^f 
and  are  the  real  and  Imaginary  parts  of  the  self 

Impedance,  and  Rnutual  ^mutual 

Imaginary  parts  of  the  mutual  Impedance.  If  the  element  Is 
operated  singly  and  not  In  an  array,  then  and 

^mutual 


Taken  alone,  pc  Is  called  the  specific  acoustic 
Impedance,  with  units  Kllogram/Meter^-Second  or  Rayl.  The 
specific  Impedances  of  air  and  water  are  415  and  1.48x10^ 
Rayls,  respectively.  It  Is  this  difference  In  specific 
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Impedance  which  accounts  for  most  of  the  variation  In 
transducer  performance.  Further,  the  transducers  are 
designed  to  be  operated  In  water,  and  the  Impedance 
mismatch  between  the  transducer  face  and  the  air  prevents 
the  transducer  ^rom  radiating  much  energy.  In  air,  the 
mechanical  dissipation  of  is  many  times  greater  than  the 
energy  radiated  Into  the  medium.* 

The  self  Impedance  Is  the  Impedance  of  the  medium  to 
the  movement  of  the  transducer  face.  It  Is  assumed  that 
the  face  Is  flush-mounted  In  an  Infinite  plane  rigid 
baffle.  The  self  Impedance  terms  are  functions  of  the  non- 
dimensional  quantity  ka ,  or  wa/c.  Thus,  the  self  Impedance 
varies  with  both  the  size  oi.  the  radiator  and  the 
frequency. 

The  mutual  Interaction  Impedance  Is  the  additional 
Impedance  on  the  transducer  face  caused  by  the  movement  of 
other  transducers  nearby.  The  mutual  impedance  varies  with 
ka  and  the  relative  spacing,  orientation,  and  velocity  of 
the  transducers.  The  total  mutual  Impedance  on  any  element 
Is  the  sum  of  contributions  from  every  other  element  In  the 
array.  Generally,  because  of  the  symmetry  of  the  array 
geometry,  more  than  one  location  In  an  array  may  have  a 
given  radiation  loading.  For  example.  In  a  five  by  five 

*  In  the  experiments  described  In  the  previous  chapter,  the 
radiation  loading  was  so  low  that  there  was  very  little 
difference  between  air  loaded  and  vacuum  (no  load) 
Immlttance  measurements. 


I 
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array,  there  are  only  four  unique  values  of  radiation 
loading,  because  of  the  sixfold  symmetry  of  a  square  array. 

More  detailed  descriptions  of  the  equations  and 
computer  programs  to  find  the  self  and  mutual  Impedance 
terms  are  given  In  Appendices  A  and  B,  respectively. 

The  calculated  radiation  impedances  can  be  accurately 
incorporated  Into  the  model  because  of  the  measurement  of 
the  transformation  factor  (|) .  Without  this  measurement,  it 
would  not  be  possible  to  predict  the  response  of  an  element 
to  the  calculated  radiation  loads.  The  alternative  is  to 
test  the  transducer  elements  with  water  loading,  which 
requires  time-consuming  bonding  techniques.  Further,  the 
effects  of  mutual  Impedance  could  only  be  seen  if  an  entire 
array  were  built  and  tested.  This  Is  not  desirable,  since 
the  goal  of  the  study  Is  to  predict  array  performance 
before  the  array  Is  built.  The  inclusion  of  radiation 
loading,  using  the  transformation  factor  ,  Is  the  key  to 
accurate  prediction  of  water-loaded  transducer  element 
performance  from  In-alr  measurements. 

4.2  The  Prediction  of  Array  Performance 

When  the  radiation  loading  for  a  given  array  location 
Is  incorporated  Into  the  lumped-parameter  equivalent 
circuit  model  of  a  transducer  element,  the  Input/output 
characteristics  of  the  circuit  can  be  used  to  predict  the 


response  of  the  element  when  It  is  used  In  that  location. 
Since  the  radiation  loading  varies  with  location,  the 
response  of  the  element  will  likewise  change  with  location. 
This  predictive  capability  can  be  used  to  select  an  optimum 
transducer  element  for  a  given  array  location.*  By 
combining  the  predicted  responses  of  every  element  chosen 
to  be  part  of  a  given  array,  it  Is  possible  to  predict  the 
response  of  the  entire  array,  before  It  Is  built. 

The  acoustic  performance  of  an  array  Is  judged  by  four 
criteria: 

1.  Transmitting  Voltage  Response  (S^)  In  dB  re  1  yPa/Volt 

2.  Transmitting  Current  Response  (S^)  in  dB  re  1  y Pa/Amp 

3.  Free-Fleld  Voltage  Sensitivity  (Mg)  In  dB  re  1  Volt/yPa 

4.  Free-Fleld  Current  Sensitivity  (M^^)  in  dB  re  1  Amp/yPa 

Sg  and  must  also  specify  the  distance  at  which  the 
measurement  Is  taken.  A  more  detailed  explanation  of  these 
terms  can  be  found  In  Albers^  or  In  the  ANSI  standard. 

These  criteria  can  be  derived  from  the  Input/output 
equations  for  the  single  transducer  element  by  the 
following  technique. 

If  an  element  k  In  an  array  Is  driven  with  some 
voltage  E,  then  the  volume  velocity  produced  by  the 
element  Is  lv|^  -  |E||Tg||^,  where  Is  the  response  of 

*  See  the  next  section. 
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the  element  as  expressed  in  Equation  5.  To  find  the 
pressure  at  a  distance  r  from  the  array,  and  hence  S^,  the 
power  radiated  by  the  element  must  be  found.  This  Is 
expressed  as  I W  |  where  (Rg)jj  is  the  real  part 

of  the  radiation  Impedance  loading  on  the  transducer 

element.  The  total  power  W  radiated  by  the  array  Is  the 
sum  of  the  W^.  The  total  transmitting  voltage  response  of 
the  array,  S^,  can  be  written  as 

Sg  -  201og[PcWD/2Trr2]  (25) 

where  P  Is  the  density  of  the  medium,  c  Is  the  sound 
velocity,  and  D  is  the  directivity  factor  of  the  array. 

The  directivity  factor  may  be  expressed  as 

D  •  4Trs/x2  (26) 

where  S  Is  the  area  of  the  array,  and  X  Is  the  wavelength. 
By  a  similar  derivation,  can  be  found  by  using  (Tj^)jj, 
the  transmitting  current  response  of  each  transducer,  as 
expressed  In  Equation  7. 

The  receiving  voltage  and  current  sensitivities, 
and  can  be  found  by  use  of  the  reciprocity  condition, 

iMgl/ISgl  •  iM^I/ISjl  •  2r/Pf  (27) 

Mg  and  are  derived  using  reciprocity  because  It  Is 
a  general  measurement  practice  to  connect  all  of  the 
elements  In  an  array  In  parallel.  Because  of  this,  the 


elements  Interact  electrically.  Therefore,  no  simple  sum 
of  Individual  element  voltages  or  current  can  correctly 
model  the  total  array  voltage  or  current. 

In  order  to  teat  the  prediction  technique,  the 
equations  for  ,  and  were  programmed  Into  the 

HP98256  calculator.  An  array  was  constructed  from  the 
elements  selected  by  the  computer  routine  discussed  In  the 
next  sections.  The  measured  frequency  response  of  the 
array  was  compared  to  that  predicted  by  the  model.  The 
results  are  shown  In  Figures  17  through  20.  The  technique 
Is  accurate  In  predicting  the  response  of  the  array,  within 
a  few  dB,  over  a  wide  frequency  range.  Some  of  the  small 
discrepancy  In  level  and  the  frequency  at  which  the  curve 
peaks  Is  probably  due  to  the  theoretical  assumptions  used 
In  calculating  the  radiation  Impedance.  Both  the  self  and 
mutual  Impedance  were  calculated  assuming  that  the  radiator 
was  In  an  Infinite  baffle,  and  the  actual  array  was  not. 

4 . 3  An  Element  Selection  Method 

The  Kendlg  scatter  diagram  method  of  element 
selection,  as  described  In  Appendix  E,  suffers  from  a 
number  of  Implementation  problems.  The  derivations  for  the 
terms  proportional  to  amplitude  and  phase  are  really  valid 
only  near  It  Is  also  assumed  In  the  derivations  that 

efficiencies  and  Qgj's  do  not  vary  significantly  from 
element  to  element.  Since  there  was  no  method  for 


FREQUENCY 


Figure  18.  Predicted  and  Measured 

Solid  lines  are  predicted  values;  dashed  lines 
are  measured  values 
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accurately  measuring  the  electromechanical  transformation 

ratio,  the  effect  of  variations  In  this  parameter  had  to  be 

Ignored.  Further,  since  this  factor  couples  In  the  effects 

of  radiation  Impedance  ,  variations  In  Impedance  can  not  be 

brought  Into  the  analysis.  It  Is  assumed  that  the  scatter 

diagram  Itself  will  show  a  Gaussian  distribution,  which  may 

not  be  the  case  for  all  groups  of  elements.  Finally,  It 

was  assumed  In  the  scatter  diagram  method  that  equal 

% 

percentage  errors  In  F^  and  (D^)  had  equal  Importance  In 
the  selection  process.  The  solutions  to  many  of  these 
problems  were  beyond  the  measurement  capabilities  available 
at  the  time.  Nonetheless,  the  technique  has  been  used  to 
construct  arrays  with  performance  better  than  arrays  built 
using  a  random  selection  method.  However,  with  the 
capability  to  predict  the  response  of  a  transducer  element 
to  any  acoustic  load,  and  over  a  range  of  frequencies,  some 
of  the  problems  with  the  Kendig  method  can  be  overcome. 

The  lumped-parameter  equivalent  circuit  model  Is  used 
as  a  flexible  tool  for  prediction  of  element  response  and 
Is  valid  over  a  broad  range  of  frequencies,  so  that  the 
element  selection  Is  not  limited  to  a  narrow  band  around 
F^ .  The  actual  amplitude  and  phase  of  element  can  be 

ij 

plotted  In  the  scatter  diagram.  Instead  of  F^  and  (D^)  • 

Of  particular  Importance  Is  the  effect  of  acoustic 
radiation  loading  on  the  response  of  the  elements.  The 
acoustic  radiation  loading  plays  an  Important  role  In 


determining  an  element's  position  In  the  scatter  diagram. 
The  relative  positions  of  elements  within  the  scatter 
diagram  change  with  loading.  As  noted  in  Appendix  B,  the 
variations  In  acoustic  radiation  loading  caused  by  mutual 
Interactions  fluctuate  over  a  very  wide  range.  The  loading 
varies  significantly  with  both  position  and  frequency. 

In  order  to  use  the  amplitude  and  phase  predictions  of 
the  model  In  a  selection  method,  the  predictions  must  be 
Incorporated  Into  some  performance  criterion.  The  Kendlg 
method  Is  essentially  a  technique  for  selecting  the  "best 
four"  elements  at  a  time  from  the  pool  of  available 
elements.  The  four  elements  are  placed  into  four  array 
locations  which  are  symmetric  about  the  center  of  the 
array.  The  most  central  array  locations  are  filled  first, 
and,  working  outward,  the  process  is  repeated  until  the 
entire  array  Is  filled. 

The  criteria  for  "best  four"  are  expressed  in  terms  of 
position  In  the  scatter  diagram.  The  criteria  are 
closeness  to  the  median  amplitude  and  phase,  proximity  of 
paired  elements  to  one  another,  and  opposition  of  pairs 
about  the  median.  Underlying  the  positional  criteria  Is  a 
derivation  of  the  beam  pattern  response  errors  due  to 

*  From  mutual  Interaction  analysis.  It  can  be  shown  that, 
for  an  unsteered  array,  four  locations  symmetric  about  the 
array  center  each  have  the  same  total  radiation  loading. 
This  will  become  Important  later  In  the  modelling  of  the 
Kend Ig  method  . 
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amplitude  and  phase  variations  among  the  four  elements. 

The  positional  criteria  were  derived  as  a  means  of 
minimizing  the  beam  pattern  errors.  The  error  derivation 
forms  a  better  mathematical  selection  scheme  since  the 
relative  Importance  of  amplitude  and  phase  deviations  tends 
to  be  distorted  In  the  scatter  diagram. 

Of  the  criteria  derived  In  the  Kendlg  method,  one  was 
chosen  for  the  present  study.  If  are  the  amplitude 
deviations,  and  are  the  phase  deviations,  then  the  total 
error  E  may  be  defined  as 

E  •  A^a^  +  A2a2  +  A3a3  +  A^a^  (28) 

This  error  criterion  was  chosen  because  of  ease  of 
Implementation,  and  because  the  other  criteria  either  did 
not  Involve  phase,  or  were  more  suited  to  selecting  pairs 
of  elements  rather  than  groups  of  four.  The  group  of  four 
elements  with  the  lowest  value  of  E  was  considered  the 
"best  f our "  . 

The  results  of  using  this  selection  technique  versus 
the  Kendlg  hand-selected  method  are  discussed  In  the  next 
sec  tlon . 

4 . 4  A  Test  Case 

The  selection  and  prediction  techniques  outlined  In 
the  previous  two  sections  were  tested  by  using  them  In  the 
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design  of  a  test  array.  The  array  was  made  up  of  the 
elements  tested  as  part  of  the  work  described  In  Chapter 
Three . 


The  design  goal  was  to  achieve  the  lowest  sldelobe 
levels  possible,  using  a  42.8  dB  shading  pattern  optimized 
by  a  Goal  Programming  technique. Naturally,  this  array 
shading  Influences  the  effects  of  mutual  Interaction 
Impedance  between  elements.  The  variation  In  element 
performance,  caused  by  different  radiation  loadings  and 
element  circuit  values,  results  In  a  distorted  array 
shading  pattern  and  higher  than  Ideal  sldelobe  levels.  It 
was  hoped  that  the  automated  selection  of  elements  would 
result  In  a  design  with  lower  sldelobes  than  a  design  using 
random  selection,  and  possibly  with  lower  sldelobes  than  a 
design  using  the  Kendlg  hand-select  method. 

The  direct  Implementation  of  a  "best  four"  approach 
to  the  selection  process  presented  a  number  of  problems. 
First  was  the  sheer  size  of  the  combinatorial  problem 
Involved  In  selecting  every  unique  set  of  four  elements 
from  a  group  of  N  elements.  For  N  greater  than  15,  the 
length  of  time  required  to  generate  and  test  the 
combinations  becomes  excessive.  Secondly,  It  was  found 
that  the  changes  In  radiation  loading  with  position  caused 
significant  shifts  In  the  amplitude  and  phase  response  of 
the  elements.  For  this  reason,  the  "four  best"  chosen  for 
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a  given  set  of  array  locations  may  not  work  well  with  the 
"best  four"  for  another  set  of  array  locations.  A  method 
was  needed  to  reduce  the  number  of  elements  under 
consideration  at  one  time  to  less  than  fifteen  and  to  try 
to  assure  that  the  elements  chosen  would  result  In  good 
performance  when  used  all  together,  and  not  just  as  sets  of 
four  . 


As  a  solution  to  the  first  problem.  It  was  decided  to 

use  the  average  element  phase  and  amplitude  as  a  reference 

point,  and  use  each  element's  mean  square  deviation  from 

this  reference  as  an  Indication  of  the  element's 

suitability.  Thus  an  element  near  the  reference  was  judged 

better  than  one  further  from  It.  To  reduce  the  number  of 

elements  considered  at  one  time,  the  element  with  the 

greatest  deviation  from  the  reference  point  was  removed  and 

a  new  average  (reference)  computed.  This  process  was 

repeated  until  less  than  fifteen  elements  remained.  Then 

the  combinatorial  subroutine  selected  the  "best  four"  from 

those  remaining.  The  four  elements  chosen  were  stored  and 

« 

removed  from  further  consideration. 

In  order  to  reduce  the  effect  of  the  shifts  In 
amplitude  and  phase,  once  the  first  four  elements  were 
chosen,  the  average  amplitude  and  phase  of  the  four  was 
stored  and  used  as  the  reference  point  for  all  subsequent 
selections.  Thus  the  process  attempted  to  find  the  "four 
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best"  elements  with  respect  to  the  average  array  amplitude 
and  phase.  This  Improved  the  uniformity  of  the  element 
response  in  the  array. 

Once  an  element  had  been  chosen  for  every  position  in 
the  array,  the  resulting  beam  pattern  was  predicted.  This 
was  done  by  finding  the  average  amplitude  and  phase  for  the 
entire  array,  and  then  expressing  the  response  of  each 
element  as  a  deviation  from  the  average.  The  ideal 
response  of  every  location  (the  shading  coefficient)  was 
multiplied  by  the  deviation  Introduced  by  the  element  used 
at  the  location.  The  result  is  a  new  array  shading 
pattern,  which  can  be  analyzed  by  standard  array  beam 
pattern  techniques.* 

The  selection  process  as  described  was  Implemented  on 
the  HP9825B  calculator.  The  program  calculated  the 
response  of  each  element  for  a  given  acoustic  radiation 
loading,  and  after  eliminating  the  elements  furthest  from 
the  reference  point,  it  selected  the  four  which  had  the 
lowest  value  of  E.  The  four  were  grouped  into  pairs  by 
matching  the  two  elements  with  the  closest  magnitude 
response.  The  four  elements  were  printed  out  and  removed 
from  further  consideration.  The  process  was  repeated, 
using  different  radiation  loadings,  until  the  array  was 
filled. 


See  Appendix  F. 
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In  order  to  gauge  the  relative  effectiveness  of  the 
selection  technique,  two  other  arrays  were  designed  for 
comparison.  For  one  array,  the  Kendlg  hand-select  scatter 
diagram  technique  was  used.  All  of  the  elements  were 
plotted  on  a  scatter  diagram,  using  and  (D^)^  as 
coordinates.  The  elements  were  selected  by  a  person 
familiar  with  the  method,  and  the  result  may  be  taken  as 
the  best  selection  possible  using  the  hand-select  method. 
Another  array  was  designed  by  taking  the  elements  selected 
by  the  automated  technique,  and  placing  them  In  the  array 
In  a  random  pattern. 

The  results  of  the  selection  routine  are  shown  In 
Figures  21  through  23.  The  lines  on  the  graph  represent: 
1)  the  Ideal  beam  pattern;  2)  the  predicted  beam  pattern 
for  the  array  designed  by  the  automated  technique;  3)  the 
predicted  beam  pattern  for  the  array  design  by  the  Kendlg 
hand-selection  method;  and  4)  the  predicted  beam  pattern 
for  the  array  of  elements  selected  at  random  from  the 
automated-design  array. 


DEGREES  OFF  ARRAY  AXIS 


Figure  21.  Comparison  of  Selection  Methods:  0°  Roll  Plane 
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22. 


Comparison  of  Selection  Methods: 


45°  Roll  Plane 


DEGREES  OFF  ARRAY  AXIS 


Figure  23. 


Comparison  of  Selection  Methods:  90°  Roll  Plane 


CHAPTER  5 


SUMMARY  AND  CONCLUSIONS 


5 . 1  Summary 

The  goal  of  this  thesis  was  to  develop  the 
experimental  and  computational  tools  necessary  to  predict 
the  performance  of  transducer  elements  In  arrays,  using  In- 
alr  measurements  of  actual  elements.  The  lumped-parameter 
equivalent  circuit  model  was  used  to  analyze  the 
characteristics  of  transducer  elements  and  predict  their 
behavior  under  varying  acoustic  loads. 

In  Chapter  1,  the  basic  problem  was  presented,  and 
expressed  In  terms  of  previous  work  which  had  attempted  to 
predict  transducer  element  behavior.  The  lumped-parameter 
equivalent  circuit  model  was  derived  In  Chapter  2,  and  the 
equations  of  that  chapter  form  the  basis  of  the  measurement 
and  prediction  techniques.  In  Chapter  3,  the  measurement 
procedure  was  developed  and  Implemented.  The  procedure  Is 
completely  automated,  and  Is  rapid,  accurate  and  repeatable 
to  within  +  2  percent.  The  self-check  routines  built  Into 
the  programming  confirm  the  accuracy  of  the  equivalent 
circuit  model  In  predicting  transducer  Immlttance  over  the 
desired  frequency  range.  The  prediction  of  transducer 
element  performance  was  presented  In  Chapter  4.  Also 
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presented  In  Chapter  4  was  a  selection  technique  for 
placing  elements  In  an  array.  The  directivity  pattern  of 
an  array  designed  by  the  selection  technique  was  compared 
to  patterns  of  arrays  designed  by  comparable  techniques, 
and  showed  5  to  lOdB  lower  sldelobe  levels,  depending  upon 
the  roll  plane.  The  prediction  of  total  array  frequency 
response  was  within  +  2dB  of  the  measured  response  of  an 
actual  array,  over  the  frequency  range  of  Interest. 

5 . 2  Conclusions  and  Future  Studies 

The  techniques  developed  In  this  thesis  can  be  used  to 
significantly  Improve  the  array  design  process,  lowering 
the  development  cost,  reducing  the  time  spent  to  design 
hardware,  and  Improving  the  quality  of  the  final  product. 
Since  array  performance  can  be  predicted  from  measurements 
of  Just  one  or  two  prototype  transducer  elements,  the  need 
to  build  and  test  entire  prototype  arrays  Is  eliminated, 
dramatically  reducing  the  time  and  cost  Involved.  With  the 
savings,  more  effort  could  be  focussed  on  testing  Improved 
designs  of  transducer  elements,  since  the  effect  of  element 
changes  on  array  performance  could  be  studied  more  easily. 
Also,  the  selection  routine  developed  has  shown  promise  as 
a  means  for  automatic  selection  of  transducer  element 
locations  within  an  array. 

The  accuracy  of  the  predictions  could  be  enhanced,  by 
elthe.  aprovlng  the  positioning  system  used  In  the 
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velocity  measurement)  or  by  attempting  to  simulate,  during 
the  Immlttance  and  velocity  measurements,  the  mounting  used 
In  array  construction. 

The  techniques  as  developed,  while  powerful  as  a 
prediction  tool,  also  form  a  solid  basis  for  further 
research . 

The  hardware  and  software  used  to  measure  the 
electromechanical  transformation  ratio  could  easily  be 
modified  to  study  the  modal  vibration  patterns  of 
transducer  elements.  Using  techniques  aleady  developed,  a 
quantitative  study  of  transducer  head  flap  and  Its 
relationship  to  transducer  resonances  could  be  made. 

The  use  of  a  technique  known  as  Non-Linear  Goal 

1  2 

Programming  or  NLGP,  as  developed  by  J.  P.  Ignlzlo, 
promises  to  further  Improve  the  element  selection 
techniques  derived  here.  By  directly  analyzing  transducer 
element  choices  In  terms  of  predicted  beam  patterns,  NLGP 
can  eliminate  most  of  the  problems  encountered  In  modelling 
the  Kendlg  method,  and  provide  a  better  final  array 
response . 


APPENDIX  A 


THE  ACOUSTICAL  IMPEDANCE  ON  A  RECTANGULAR  PISTON 


The  expression  for  the  acoustical  radiation  Impedance 

1  3 

Is  derived  In  Morse  &  Ingard.^-*  In  section  7.4,  which  deals 
with  radiation  from  a  plane  surface,  the  Green's  function 
Is  solved  to  find  the  reaction  force  on  a  rectangular 
piston  with  sides  of  length  a  and  b.  The  expression  Is 
simplified  by  comparison  of  series,  and  for  the  case  of  a 
square  piston  (a»b),  the  solution  reduces  to* 

Eq.  7,4.44  =  pcu^ab [ 0 p( ka)  +  jX^Cka)]  (29) 

where  F  and  u  are  the  force  on  and  velocity  of  the  piston 

OJ  U) 

face,  and  6^  and  X^  are  given  by 

Qq  »  1  -  2J^(ka)/ka 

Eq.  7.4.31  (30) 

IT  /  2 

Xo  *  Mj^(ka)  -  4^  8in(  (  ka)  cos(  q)  )  sln^  (  q)  dq 

Is  the  first  order  Bessel  function,  and  Is  the  first 
order  Struve  function. 

The  computation  of  on  a  large  computer  Is 

*  This  equation  has  been  rewritten  to  conform  to  the 
exponential  time  dependence  convention  (e^-^  ^)  used  In  this 
study. 
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straightforward,  since  the  answer  can  be  obtained  through 
standard  scientific  subroutine  packages.  The  Struve 
function  requires  numerical  Integration  which  was  done 
using  a  three  point  Gaussian  quadrature  subroutine .  The 
accuracy  of  the  routine  depends  upon  the  number  of 
subintervals  taken  within  the  limits  of  Integration.  The 
Struve  function  converges  fairly  rapidly,  and  seven  places 
of  accuracy  were  achieved  with  little  computation  time. 

These  expressions  for  the  radiation  loading  were 
Implemented  In  the  subroutine  RADIMP.  The  subroutine 
returns  the  values  of  9^  and  for  a  square  radiator  used 
In  water,  given  the  length  of  one  side  and  the  frequency  of 
operation.  The  accuracy  of  the  subroutine  was  checked  by 
comparing  results  with  Morse  and  Ingard's  Table  IX. 

The  subroutine  was  used  within  a  larger  program  which 
found  both  the  self  and  mutual  Impedance  for  an  array  of 
square  radiators.  The  program  Is  listed  In  Appendix  G. 


APPENDIX  B 


THE  MUTUAL  IMPEDANCE  BETWEEN  RECTANGULAR  PISTONS 


The  expression  for  Che  acoustical  mutual  Interaction 
Impedance  between  two  rectangular  pistons  is  derived  by 
Arase.^^  Figure  24  gives  the  physical  layout  of  two 
Identical  radiators,  and  the  dimensions  used  in  Che 
derivation.  If  we  let 


Eq.  9  ka  ■  A  kb  «  B  kg  •  G  kh  ■  H  (31) 

then  the  mutual  radiation  resistance  becomes 


Eq.  10  (Rjjj/Pcab)  •  (  1 /2TrAB )  {  4C  ( H  ,  G)  -  2C(H-A,G) 

-  2C(H,G-B)  -  2C(H+A,G)  -  C(H-A,G+B) 

+  C(H+A,G-B)  +  C(H+A,G+B)}  (32) 


where 


Eq.  11  C(x,y)  ■  |x| 


(y^+x^)*^ 


[(r^  -  y^)^(co8(r)/r)]dr 


+  lyl 


|y| 

(y^+x^)*^ 


[(r^  -  x^)^ (co8(r)/r) ]dr 


|y| 


-  cos 


(x^  +  y2)^  -  (x^  +  y2)^sln(x2  + 


(33) 


The  expression  for  the  mutual  radiation  reactance 
divided  by  pcab  has  the  same  form  as  Eq.  32,  but  the 
function  C(x,y)  Is  replaced  by  S(x,y)  In  each  case,  where 
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Eq .  12 


S(x,y) 


(y^+x^)*^ 
1x1  I  [(r2 

\y\ 

(y^x^)^ 


-  lyl 


[(r‘ 


|y| 


x^)^(sln(r)/r) ]dr 


x^)**(sln(r)/r)  ]  dr 


-  (x' 


y2)^cos{x^ 


y2)^+  sln(x2 


2^h 


+  y^) 


(34) 


The  computation  of  and  X^^  Involved  the  numerical 
Integration  of  Eqs.  33  and  34,  using  the  same  three  point 
Gaussian  quadrature  routine  discussed  In  Appendix  A.  The 
listing  of  the  program  to  find  the  total  loading  on  an 
8X8  array  of  elements  Is  given  In  Appendix  G. 


The  total  mutual  Interaction  Impedance  on  an  element 
In  an  array  Is  the  weighted  sum  of  the  contributions  from 
every  other  element  In  the  array.  If  Is  the  Impedance 

on  the  element  due  to  the  movement  of  the  element, 

then  the  total  loading  on  the  element  In  an  array  of  M 

elements  Is 

M 

"  Z  (Uk/Uj)ZjJ.k  (35) 

k-1 


75 


where  Uj^  and  Uj  are  the  velocities  of  the  and 

elements  respectively,  and  +  jXjjj.  If  we  label  the 

self  Impedance  on  the  element,  as  found  In  the  previous 
Appendix,  as  ,  then  the  total  loading  on  the  element 

1  s 

Zjotal  -  ZJ  +  I  (Uk/Uj)zJ*>^  (36) 

k=l 

From  this  It  Is  evident  that  the  loading  on  an  element 
wll  vary  with:  1)  position  In  the  array,  and  2)  the  array 
shading  distribution,  which  affects  the  relative  velocities 
of  the  elements.  For  the  case  of  an  8  X  8  array,  as  shown 
In  Figure  25,  there  are  ten  unique  values  of  radiation 
loading,  due  to  the  symmetry  of  the  array.  The  letters  A 
through  J  denote  those  positions. 

To  Illustrate  the  first  point,  the  variation  of 
loading  with  position  and  frequency  Is  shown  for  locations 
A,  D,  and  H  In  Figures  26  and  27.  Figure  26  Is  for  the 
real  part  of  the  loading;  Figure  27  is  for  the  imaginary 
part  of  the  loading.  The  solid  line  Is  the  loading  on  a 
single  element  and  the  line  on  the  frequency  axis  denotes 
the  frequency  for  which  k  ■  it  .  It  Is  easily  seen  that  the 
loading  varies  significantly  with  position  and  frequency, 
especially  below  k  ■  if  . 


I 

I 


maglnary  Radiation  Loading  for 
ay  Positions 
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Table  1  demonscraces  the  variation  in  loading  caused 
by  changing  the  array  shading.  The  values  in  the  table 
have  been  normalized  by  the  loading  on  a  single  element. 
The  frequency  used  for  this  example  was  chosen  such  that 
k  «  IT. 


r 


TABLE  1 


Comparison  of  Radiation  Loading  by  Position 
For  Various  Shading  Patterns 


POSITION 

LOADING* 

UNSHADED 

SHADING  1 

SHADING  2 

A 

REAL 

1.041 

1.082 

1.082 

IMAG 

0.063 

0.126 

0.131 

D 

REAL 

0.975 

1.077 

1.095 

IMAG 

0.017 

0.105 

0.106 

E 

REAL 

1.232 

1.130 

1.122 

IMAG 

-0.033. 

0.097 

0.093 

F 

REAL 

1.176 

1.157 

1.161 

IMAG 

0.282 

1 . 372 

0.055 

H 

REAL 

0.970 

1.046 

1.124 

IMAG 

0.370 

0.057 

0.057 

POSITION  SHADING  COEFFICIENTS  FOR  SHADING 

1 

A:  1.0 

B:  0.7606 

C:  0.5844 

D:  0.4672 

E: 

0.3561 

F:  0.1882 

G:  0.2196 

H:  0.1186 

I:  0.0102 

J: 

0.0015 

POSITION  SHADING  COEFFICIENTS  FOR  SHADING 

2 

A:  1.0 

B:  0.7282 

C:  0.5263 

D;  0.3815 

E: 

0.2632 

F:  0.1095 

G:  0.1189 

H:  0.0612 

I:  0.0131 

J: 

0.0014 

’^NOTE:  The  radiation  loadings  have  been  scaled  as  explained 
in  the  text> 


APPENDIX  C 


USING  THE  HP3750  TO  MEASURE  IMMITTANCE 


The  HP3570  Network  Analyzer  Is  designed  to  measure  the 
transfer  function  response  of  two-port  networks.  However, 
with  the  use  of  special  Input  terminations  and  simple 
formulae,  the  analyzer  can  be  used  to  find  the  Input 
Immlttance  of  one-port  devices. 

The  standard  mode  of  operation  of  the  analyzer  Is 
shown  In  Figure  28.  The  analyzer  supplies  a  voltage  V  to 
two  output  terminals  and  which  are  connected  to 

the  Inputs  of  the  networks  under  test.  Generally,  one 
network,  fed  by  A^^^ ,  is  used  as  a  reference  or  standard, 
and  the  network  fed  by  Is  to  be  compared  with  the 

reference.  The  outputs  of  the  reference  and  test  networks 
are  connected  Into  Aj^^j  and  respectively.  The  analyzer 

measures  A^^  and  B^^,  and  finds  the  ratio  E  ~  I  A^jj  I  / 1  B^,,  I 
and  the  phase  angle  P  ■  ArgCA^^j)  -  Arg(Bj^^).  The  resistors 
marked  R^^^^  are  the  Internal  source  Impedance  of  the 
driving  voltage  sources,  and  the  resistors  marked  are 

external  resistors  used  for  Impedance  matching. 

To  use  the  analyzer  for  immlttance  measurements,  the 
set-up  shown  In  Figure  29  Is  used.  The  A  channel  Is 
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shorted,  except  for  the  matching  resistor  Rg^t’  ® 

channel  is  connected  to  the  unknown  Impedance  Z,  again  with 
the  matching  resistor. 


The  equivalent  circuits  for  the  equipment  set-up  in 
Figure  29  are  shown  in  Figure  30.  Solving  for  and 


^In 

Bin  ■  V{Rext/(Rint  +  ^ext  +  Z)} 


(37) 


Since  the  unknown  Impedance  can  be  both  resistive  and 
reactive,  let  Z  »  A  +  jX.  Substituting  into  the  above  and 
solving  for  E  and  P  yields 


(38) 

P  .  tan-Mx/(Ri„^  +  Rg^,  +  R)} 

These  two  equations  can  be  solved  for  R  and  X,  which  yields 


+  Rext>^Bco8(P)  -  1] 

^  -  (Bint  +  Rext>B^in(P) 


(39) 


For  the  experiments  described  in  Sections  3.1  and  3.2,  the 
resistors  Rj^t  Bg^t  equal  to  75  Ohms.  Therefore 


R  ■  150[Ecos(P)  -1]  X  •  150E8in(P) 


(40) 
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Solving  for  the  admittance  yields 

G  -  {co8(P)E  -  l)/{l50tE^  -  2co8(P)E  +  1]*^} 

B  •  8ln(P)E/{150[E2  -  2co8(P)E  +  I]'®} 


(41) 


These  last  two  sets  of  equations  were  programmed  Into 
the  HP9825B  calculator  to  find  the  Immlttance  of  the 
transducer  elements. 


APPENDIX  D 


THE  INFLUENCE  OF  Rq  ON  THE  EQUIVALENT  CIRCUIT 

In  order  to  assess  the  Influence  of  R^  on  the 
equivalent  circuit,  the  following  method  was  used  to 
estimate  R^  from  the  measurements  described  In  Section  3.1. 

From  Figure  13,  It  was  assumed  that  the  difference 
between  the  measured  and  predicted  Impedance  peaks  at  ^ 
was  due  to  neglecting  R^  In  the  equivalent  circuit. 

Figure  31  shows  how  R^  would  affect  the  equivalent  circuit. 
If  R2  Is  the  predicted  resistance  of  the  equivalent  circuit 
evaluated  at  then  R2  would  be  lower  than  R2  >  depending 

upon  the  value  of  R^ .  If  R^  were  set  to  the  measured 
resistance  of  the  transducer  element  at  then  the 

circuit  can  be  solved  to  find  R^ .  Solving  for  R^  yields 

RoI‘^2  "  ‘ 

Rq  varies  Inversely  with  frequency,  and  the  r.  ^nce 
Is  more  properly  expressed  as  Rq  ■  (Rq/**’).  where  Rq  1  e 
theoretical  DC  resistance.  Multiplying  the  value  of  Rq 
found  using  Eq.  42  by  gives  an  estimate  of  Rq  . 

For  the  transducer  elements  measured  In  this  study, 
the  values  of  Rq  found  in  this  manner  ranged  between 
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Rz'  = 


Figure  31.  Effect  of  Adding  to  the  Equivalent  Circuit 


1.5  and  3.0  X  10^^  Ohms.  When  this  was  put  into  the  proper 
equivalent  circuit  expressions,  the  difference  In  transmit 
and  receive  levels  did  not  vary  more  than  .1  dB. 

Therefore,  It  was  concluded  that,  for  the  transducer 
elements  tested,  the  ommlsslon  of  did  not  significantly 


APPENDIX  E 


THE  KENDIG  SCATTER  DIAGRAM  SELECTION  METHOD 


In  his  report,  (Ref.  [1]),  Kendlg  describes  a  method 
for  reducing  the  detrimental  effects  of  non-ideal  elements 
on  array  performance.  He  developed  a  selection  method  for 
the  placement  of  elements  within  an  array,  based  upon  their 
relative  amplitude  and  phase  response.  The  first  part  of 
his  report  Is  a  derivation  of  the  parameters  needed  for 
prediction  of  element  response;  the  second  part  Is  a 
description  of  a  selection  scheme  for  choosing  element 
positions  In  an  array  to  reduce  minor  lobe  levels. 

Kendlg  shows  that,  near  ,  the  free  field  voltage 
sensitivity  M^  Is  approximately  proportional  to  (nD^)  , 

where  n  Is  the  efficiency  of  the  element,  and 

■  *^effQz/^<C  +  C^)(l  -  k2^f)  (43) 
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these  parameters  did  not  vary  significantly  from  element  to 
element.  In  another  report, Kendlg  had  shown  that  phase 
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errors  in  hydrophone  elements  were  proportional  to 
assuming  that  the  positioning  of  the  elements  In  the  array 
could  be  held  to  close  tolerance.  Thus,  by  comparing  the 

h 

values  of  and  (D^)  »  a  relative  measure  of  transducer 

element  amplitude  and  phase  could  be  made  from  In-alr 

measurement  procedures.  The  comparison  was  made 

graphically  In  the  form  of  a  scatter  diagram,  such  as 

Figure  32.  Each  point  represents  a  transducer  element;  the 

X  represents  the  median  combination  of  F_  and  (D_)  . 

z  z 

Kendlg  assumed  that  the  distribution  of  amplitude  and  phase 
errors  would  follow  a  Gaussian  distribution.  The  closeness 
of  an  element  to  the  median  gives  an  Indication  of  the 
relative  magnltue  of  amplitude  and  phase  errors. 

Next,  Kendlg  described  the  effect  of  amplitude  and 
phase  errors  on  array  performance.  Figure  33  represents 
two  pairs  of  array  elements  which  are  symmetrically  spaced 
about  the  array  center.  If  Is  the  amplitude,  and  a^^  Is 
the  phase  error  of  an  element,  then  the  contribution  of 
these  four  elements  to  the  horizontal  array  pattern  may  be 
expressed  as 

P(real)  "  [(A2+A2)  +  (A3+A^) ] cos( nu)  - 

[  (A]^a^-A3a3)  -  (A2a2-A4a^  )  ]  sln(  nu)  (44) 

^(Imag)  *  [(A1-A3)  +  (A2“A4 ) ] sln( nu)  + 

[(A^a^+A3a3)  +  (A2a2''-A4a^  )  ]  cos(  nu) 


(45) 
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where  n  Is  an  Integer,  u  Is  equal  to  da)sln(6) /2c ,  d  Is  the 
element  spacing,  and  6  is  the  Incident  angle.  In  the  Ideal 
case,  all  of  the  A^s  are  equal,  and  all  of  the  a^^s  are 
zero.  The  terms  in  the  above  equations  are  grouped  to  show 
the  offsetting  contributions  from  specific  pairs  within  the 
group  of  four.  It  is  desired  that  both  the  second  term  of 
Eq .  44  and  all  of  Eq .  45  be  equal  to  zero.  Further,  the 
first  term  of  Eq .  44  should  be  equal  to  four  times  the 
median  element  amplitude. 

Applying  the  two  equations  above  to  the  scatter 
diagram,  Kendlg  proposed  the  following  selection  process: 
pick  two  elements  which  are  close  to  one  another  in  the 
scatter  diagram,  such  as  A  and  B,  and  place  them  in 
positions  1  and  3.  By  doing  this,  the  (A^a]^-A3a3)  and 
(Aj^-A3)  terms  are  made  small.  Elements  C  and  D  are  then 
chosen  for  positions  2  and  4.  These  two  elements  are  also 
close  together  in  the  scatter  diagram,  which  makes 
(A2a2~A^a^}  and  (A2~A^)  small.  Further,  since  C  and  D  are 
nearly  opposite  in  amplitude  and  phase  error  from  A  and  B, 
the  second  term  of  Eq .  45  will  be  small,  and  the  first  term 
of  Eq .  44  will  be  nearly  equal  to  four  times  the  median 
element  amplitude. 


APPENDIX  F 

ARRAY  BEAM  PATTERN  RESPONSE 


The  beam  pattern  response  of  an  array  Is  the  sum  of 
the  contributions  from  all  of  the  elements  within  the 
array.  The  following  derivation  Is  for  a  planar  array  of 
point  elements,  but  Is  can  easily  be  extended  to  find  the 
response  of  an  array  of  square  radiators. 

Consider  a  point  element  m  In  a  plane  array  located  at 
position  shown  In  Figure  34.  For  a  wave  coming 

from  the  direction,  the  contribution  of  the  m^*’ 

element  to  the  total  array  signal  Is 


Qm 


A_e 

m 


j(a)t  +  u) 


"  Aj,eJ“^[  <cos(u)  +  jsln(u)] 

where 


(46) 


u  ■  kX|||Sln6cos(^  +  kyu|Sln(|>sln(})  +  3^ 

and  A|||  Is  the  element  amplitude  shading  and  a^j  Is  the 
element  phasing.  The  e^^’^  time  dependence  may  be  dropped, 
since  It  Is  common  to  all  the  elements,  and  therefore  does 
not  contribute  any  Information  to  the  directional  response. 
The  total  response  of  the  array  of  M  elements  Is  simply  the 


sum  of  the  contributions  of  each  element 


NORMAL  TO  INCIDENT 
WAVE 


V„) 
tn  m 
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Qtot  “  Z  Am{  co8[kx„8ineco84)  +  kyaBin08la(|)  +  a^] 

m=i 

+  j8ln[  kX|n8ln6co8({>  +  ky^Bl  a6  8  ln(|>  +  a^,]  }  (47) 

The  aubrouClne  PRESS  waa  uaed  to  Implement  Eq.  47, 
returning  the  array  amplitude  and  phase  for  any  direction. 
For  the  beam  pattern  predictions  described  In  Section  4.3, 
each  Ideal  and  a^  was  modified  to  Include  the  error 
Introduced  by  the  m^^  element. 

For  an  array  of  square  elements  of  side  C,  the 
directivity  of  the  elements  must  be  Included  In  Eq.  47. 

This  Is  done  by  modifying  Aq,: 

Am  ■  Ajii  [  3  ln(  .  5kC  I  sln6  slncj)  I  )/ (  .  5kC  I  s  In©  co8(J)  I  )  ] 

X  [  s ln(  .  5kC  I  s ln9co8(J)  I  ) /  (  .  5kC  1  8 In©  co8<J>  I  )  ]  }  (48) 

A  listing  of  the  subroutine  PRESS  follows. 
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SUBROUTINE  PRESS(THETAP , PHIP , AMPL, PHASE) 

Q  ********************************************************* 

C 

C  SUBROUTINE  TO  FIND  ARRAY  RESPONSE 

C 

C  RETURNS  AMPLITUDE  AND  PHASE  (DEG)  RESPONSE  OF  AN  ARRAY. 
C 

C  INPUTS  ARE  BEAM  ANGLE  THETA  AND  ROLL  PLANE  PHI ( DEGRE ES) . 

C  ELEMENT  LOCATIONS  ARE  IN  ARRAYS  X  AND  Y; 

C  ELEMENT  SHADINGS  ARE  IN  THE  COMPLEX  ARRAY  ELEMNT, 

C  STORED  IN  THE  THE  FORM  OF  MAGNITUDE  AND  PHASE  (DEG). 

C  AK  IS  THE  WAVENUMBER  (K)  AND  M  IS  THE  NUMBER  OF  ELEMENTS. 
C 

Q  ********************************************************* 

COMMON  X( 101) ,Y(101) .ELEMNT (101) ,AK,M 
REAL*4  ANGLE(101),A(101) 

REAL*8  U ,SUMR,SUMI,DC0S,DSIN,DATAN2,DAMPL 
COMPLEX  ELEMNT 

C  CONSTANTS  TO  CONVERT  DEGREES  TO  RADIANS  AND  VICE  VERSA 
DATA  RAD2D/57. 2958/, D2RAD/1. 74532925  E-02/ 

C  INITIALIZE  SUMS  AND  SET  ANGLE  PARAMENTER 
SUMR  -  SUMI  -0.0 

UCOS-SIN(THETAP*D2RAD)*COS(PHIP*D2RAD) 

USIN«SIN(THETAP*D2RAD)*SIN(PHIP*D2RAD) 

C  SUM  OVER  ELEMENTS 
DO  10  J-1,M 
A(  J)-REAL(ELEMNT(J)) 

ANGLE( J)-AIMAG(ELEMNT( J) ) 

U-DBLE(AK*X(J)*UCOS+AK*Y( J)*USIN+ANGLE(J)*D2RAD) 
DAMPL-DBLE(A( J) ) 

SUMR-SUMR+DCOS(U)*DAMPL 
10  SUMI-SUMI+DSIN(U)*DAMPL 

C  TAKE  MAGNITUDE  OF  AMPLITUDE  AND  PUT  PHASE  IN  DEGREES 
AMPL-SNGL(SUMR*SUMR+SDMI*SUMI) 

AMPL-SqRT(AMPL) 

PHASE-SNGL(DATAN2 (SUMI, SUMR) )*RAD2D 

RETURN 

END 


APPENDIX  G 


LISTING  OF  THE  RADIATION  IMPEDANCE  PROGRAM 


The  following  program  was  used  to  find  the  total 

radiation  loading  on  the  array  of  elements. 

C  PROGRAM  TO  FIND  THE  TOTAL  LOADING  ON  THE 
C  UNIQUE  ARRAY  POSITIONS  IN  AN  8X8  ARRAY. 

C  THE  RADIATION  LOADING  IS  THE  COMBINATION 
C  OF  THE  SELF  IMPEDANCE,  FOUND  IS  SUBROUTINE 

C  RADIMP,  AND  THE  MUTUAL  INTERACTION  IMPEDANCE, 

C  FOUND  IN  SUBROUTINE  MUTRAD. 

C  THE  MUTUAL  IMPEDANCE  COEFFICIENTS  ARE  GENERATED 

C  BY  SUBROUTINE  FINDMT. 

C  THE  OUTPUT  VARIES  DEPENDING  UPON  THE  SHADING  PATTERN. 

C  THE  INPUTS  TO  THE  PROGRAM  ARE  THE  LENGTH  OF  THE 

C  SIDE  OF  AN  (SQUARE)  ELEMENT,  AND  THE  STARTING 
C  FREQUENCY,  THE  NUMBER  OF  FREQUENCY  STEPS  AND 
C  THE  SIZE  OF  THE  FREQUENCY  STEPS. 

IMPLICIT  REAL*8  (E,M,X,Y,V) 

REAL*8  ELEMNT(52,3) ,MUTARA(35,4) ,FREQ(21) , LENGTH 
REAL*8  DKA,RM,XM,RMUTOT,XMUTOT,RTOT,XTOT, VELRAT 
REAL*8  UNIPOS(8,3),LOAD(21,9,3) 

DIMENSION  ARAPOS(8) ,M1(4),M2(4) 

C  SETTING  UP  ARRAY 

COMMON/MUTUAL/MUTARA 

DATA  Ml/1 .  ,-l 

DATA  M2/1 . ,1 .  ,-l  .  ,-l  ./ 

DATA(UNIPOS( J,l) ,J-1 ,8)/0. 5D0, 2*1 .5DO, 2*2 .5D0, 
3.5D0,2.5D0,3.5D0/ 

DATA(UNIPOS(J,2) ,J-l,8)/2*. 5D+0, 1 .5D+0, .5D+0,1.5D+0 
1 , . 5D+0 , 2 . 5D+0 , 1 . 5D+0/ 

DATA(UNIPOS( J,3) , J-1 , 8) / 8* 1 . OD+0/ 

DATA  (ARAPOS(J) ,J-1 ,8)/lHA,lHB  ,1HC  ,1HD,1HE, 1HF,1HG, IHH/ 
K-0 

DO  27  J-1, 8 
DO  25  N-1,4 
K-K+1 

ELEMNT(K,l)-Ml(N)*DNIPOS(J,l) 

ELEMNT(K,2)-M2(N)*ONIPOS( J,2) 

ELEMNT(K,3)-UNIPOS(J,3) 

IF(UNIPOS(J,l).EQ.UNIPOS(J,2))  GO  TO  25 
K-K+1 

ELEMNT(K,l)-Ml(N)*UNIPOS( J,2) 

ELEMNT(K,2)-M2<N)*ONIPOS( J,l) 

ELEMNT(K,3)-UNIPOS(J,3) 

25  CONTINUE 
27  CONTINUE 


o  o 
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ENTER  LENGTH,  STARTING  FREQ.,  FREQ.  STEP  SIZE,  AND 
NUMBER  OF  STEPS 
READ(5,110)  LENGTH 
READ(5,lll)STRTFQ,STEPFQ,NUMSTP 
DO  30  N  -  1,NUMSTP 

FREQ(N)  -  STRTFQ  +  STEPFQ*FLOAT( N-1 ) 

C  FIND  SELF  IMPEDANCE,  AND  FILL  MUTUAL  IMPEDANCE  ARRAY 
CALL  RADIMP(FREQ(N) .LENGTH, RM,XM) 

CALL  FINDMT(FREQ(N) .LENGTH) 

C  SELF  IMPEDANCE  IS  IN  THE  FIRST  OUTPUT  ARRAY  POSITION 
LOAD(N,l,l)  -  (.15D+04)*(.1026D+04)*RM/(LENGTH**2.) 
L0AD(N,1,2)  -  (.15D+04)*(.1026D+04)*XM/(LENGTH**2.) 

C  SUM  MUTUAL  IMPEDANCES  OVER  ALL  ELEMENTS 
DO  20  L-1,8 
RMUTOT  -  O.OD+00 
XMUTOT  -  O.OD+00 
DO  15  K-1,52 

1F(ELEMNT(K,1).EQ.UNIP0S(L, 1). AND. 
1ELEMNT(K,2).EQ.UNIP0S(L,2))  GO  TO  15 
XDIFF  -  DABS(UNIPOS(L,l)-ELEMNT(K,l)) 

YDIFF  -  DABS(UNIPOS(L,2)-ELEMNT(K,2) ) 

VELRAT  -  ELEMNT(K,3)/UNIPOS(L, 3) 

CALL  MUTRAD (XDIFF, YD IFF ,RMUT ,XMUT) 

RMUTOT  -  RMUTOT  +  VELRAT*RMUT 
XMUTOT  -  XMUTOT  +  VELRAT*XMUT 
15  CONTINUE 

C  ADD  SELF  AND  MUTUAL  IMPEDANCES 
RTOT  -  RMUTOT  +  RM 
XTOT  -  XMUTOT  +  XM 

LOAD(N,L+l,l)  -  (.15D+04)*(.1026D+04)*RT0T/(LENGTH**2.) 
LOAD(N,L+l,2)  -  (.15D+04)*(.1026D+04)*XTOT/(LENGTH**2.) 
CONTINUE 
CONTINUE 
WRITE(6,113) 

WRITE(6 , 114) 

WRITE(6 , 115) 

WRITE (6, 116) (ARAPOS(J) , J-1 ,8) 

DO  40  N-1, 21 

WRITE(6, 120)1 DINT(FREQ(N) ) , ( LOAD( N, J , 1 ) , J-1 , 9) 

WRITE(6 ,121) (LOAD(N, J,2) , J-1 , 9) 

F0RMAT(F10.6) 

FORMAT(2F10.2,I2) 

F0RMAT(5X,'  ',///) 

F0RMAT(51X,'P  OSITION  IN  THE  ARRA  Y') 
FORMAT(13X,' SINGLE') 

FORMAT ( IX, 'FREQUENCY' , 3X, ' ELEMENT' , 8X ,A1 , 7( 1 IX, A1 ) , /) 
F0RMAT(1X,'FREQUENCY',18X,A1,7(11X,A1 ),////) 
FORMAT(3X,I5, 1X,9( 3X,D9. 3) ) 

F0RMAT(9X, 9( 3X,D9. 3) , /) 

CONTINUE 

WRITE(6,117)(ARAP0S(J) ,J-1,8) 

STOP 
END 

C  SUBROUTINE  TO  FIND  THE  MUTUAL  INTERACTION  IMPEDANCE 
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C  COEFFICIENT  BETWEEN  THE  REFERENCE  ELEMENT  AND  THE 
C  ELEMENT  AT  THE  POSITION  PASSED  TO  THE  SUBROUTINE. 
SUBROUTINE  MUTRAD( XDIFF , YD  IFF , RMUT ,XMUT) 

IMPLICIT  REAL*8  (E,M,X,Y,V) 

REAL*8  MUTARA(35,4) 

COMMON /MUT  UAL /MUTAR  A 
DO  10  K-1,35 

IF ( (DABS ( MUT ARA(K,1) -XDIFF) .LE  .0  .  lD-0  3  .  AND  . 
1DABS(MUTARA(K,2)-YDIFF) . LE . 0 . lD-3 ) . OR  . 
2(DABS(MUTARA(K,1)-YDIFF) .LE.0.1D-3.AND. 
3DABS(MUTARA(1C,2)-XDIFF)  .LE.O.lD-3)  )  GO  TO  20 
10  CONTINUE 

WRITE(6,100)  XDIFF, YDIFF 

100  F0RMAT(1X,'N0  MATCH  FOR  ',F5.2,'  ,  ',F5.2,'  CASE.',/) 
RMUT  -  XMUT  -  0.0 
RETURN 

20  RMUT  -  MUTARA(K,3) 

XMUT  -  MUTARA(K,4) 

RETURN 

END 

C  SUBROUTINE  TO  FIND  THE  RADIATION  IMPEDANCE  PARAMETERS 
C  THETA  NOUGHT  AND  XI  NOUGHT  AS  USED  IN  EQUATION  7.4.44 
C  OF  MORSE  AND  INGARD'S  BOOK,  'THEORETICAL  ACOUSTICS' 
SUBROUTINE  RADIMP( FREQ , LENGTH, THNOT , XI NOT) 

REAL*8  GAMMA, THNOT, XIN0T,MMBSJ1 ,PI,ANS,FREQ  , LENGTH 
DATA  SSPEED/ .15D+04/ 

PI  •  4.0*DATAN(1 .ODOO) 

GAMMA  -  2.  *  PI  *  FREQ  *  LENGTH  /  SSPEED 
CALL  INTGRl (O.ODOO ,Pl/2. , 3 , 20 , GAMMA , ANS ) 

XINOT  -  4.0  *  ANS/PI 

THNOT  -  l-(2./GAMMA)*MMBSJl(GAMMA,IER) 

RETURN 
END 

SUBROUTINE  TO  FIND  THE  MUTUAL  INTERACTION  COEFFICIENTS 
FOR  ANY  COMBINATION  OF  POSITIONS. 

THE  OUTPUT  IS  PUT  IN  THE  ARRAY  MUTARA, 

WHICH  CONTAINS  THE  X-OFFSET,  Y-OFFSET,  REAL,  AND 
IMAGINARY  INTERACTION  COEFFICIENTS. 

SUBROUTINE  F INDMT ( FREQ , LE NGTH) 

REAL*8  CA,CB,CG,CH,C , A, B,G,H,MUTARA( 35,4) , SSPEED 
REAL*8  FREQ, LENGTH 
COMMON/MUTUAL/MUTARA 
DATA  SSPEED/ .15D+04/ 

C  -  8.  *  DATANd.O  DO)*FREQ/SSPEED 
A  -  B  -  1  .0 
H  -  G  -  0.0 
L  -  0 

DO  6  1-2,8 
DO  5  K  -  1,1 
L  -  L  +  1 

G  -  MUTARA(L,1)  -  DFLOAT(I-l) 

H  -  MUTARA(L,2)  -  DFLOAT(K-l) 

CA  -  C  *  A*LENGTH 
CB  -  C  *  B*LENGTH 
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CG  -  C  *  G*LENGTH 
CH  -  C  *  H*LENGTH 

CALL  MUTIMP(CA,CB,CG,CH,MUTARA(L, 3) ,MUTARA(L, 4) ) 

5  CONTINUE 

6  CONTINUE 
RETURN 
END 

C  SUBROUTINE  TO  FIND  THE  REAL  AND  REACTIVE  PARTS  OF  THE 
C  MUTUAL  IMPEDANCE  RM  AND  XM ,  GIVEN  THE  ALPHA , BET A , G ,  AND  H 
C  AS  DESCRIBED  IN  ARASE. 

SUBROUTINE  MUT IMP ( CKA , CKB , CKG , CKH ,R , X ) 

REAL*8  CKA,CKB,CKH,CKG,R,X,CONSTl(9),CONST2(9 ), 
1C0NST3(9) ,TOTC,TOTS,PI ,C,S 
DATA  C0NSTl/4.0,4*-2.0,4*l .0/ 

DATA  CONST 2/0. 0,-1 .0,0.0,1.0,0.0,2*-1.0,2*1.0/ 

DATA  CONST3/ 2*0. 0,-1 .0,0. 0,1. 0,-1. 0,1. 0,-1. 0,1.0/ 

PI  -  4.*DATAN(1.0  DO) 

TOTS  -TOTC  -  0.0  DO 
DO  1  K-1,9 

CALL  CANDS(CKH+CONST2(K)*CKA,CKG+CONST3(K)*CKB,C,S) 

TOTC  -  TOTC  +  C*C0NST1(K) 

1  TOTS  -  TOTS  +  S*C0NST1(K) 

R  «  T0TC/(2 .*PI*CKA*CKB) 

X  •  T0TS/(2.*PI*CKA*CKB) 

RETURN 

END 

C  SUBROUTINE  TO  FIND  BOTH  C(X,Y)  AND  S(X,Y)  AT  THE  SAME  TIME 
SUBROUTINE  CANDS(A,B ,C ,S) 

REAL*8  A,B,C,S,QUANTP,ABSA,ABSB ,DC0S,DSIN,0UT(2) , 
1PLMULT(2) ,0UT1 ,0UT2 
DATA  PLMULT/1.0  ,-1.0/ 

QUANTP  -  DSQRT(A  *  A  +  B  *  B) 

ABSA  -  DABS(A) 

ABSB  -  DABS(B) 

N-15 

DO  1  K-1,2 

CALL  INTGR1(ABSB,QUANTP,N,K,B,0UT1) 

CALL  INTGRl(ABSA,qUANTP,N,K,A,0UT2) 

OUT(K)  -  PLMULT(K)*(ABSA*0UT1  +  ABSB*0UT2) 

1  CONTINUE 

C  -  OUT(l)  -  DCOS(QUANTP)  -  QUANTP*DSIN(QUANTP) 

S  -  0UT(2)  -  qUANTP*DCOS(QUANTP)  +  DSIN(qUANTP) 

RETURN 

END 

C  INTEGRAND  FROM  MORSE  AND  INGARD'S  Eq  7.4.44  AND  ARASE  PAPER 
DOUBLE  PRECISION  FUNCTION  CKERNL ( q , N , PAR AM) 

REAL*8  q,PARAM,qUANTM, check 
CHECK  -  Q  *  q  -  PARAM  *  PARAM 
IF  (CHECK. GE. 0.0  DO)  GO  TO  10 
CHECK  -  0.0  DO 
10  qUANTM  -  DSqRT( CHECK) 

IF  (FLOAT(N)  -  2.)  1,2,3 
C  INTEGRAND  FOR  C(X,Y) 

1  CKERNL  -  DCOS(q)  *  qUANTM  /  q 
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RETURN 

C  INTEGRAND  FOR  S(X,Y) 

2  CKERNL  -  DSIN(Q)  *  QUANTM  /  Q 
RETURN 

C  INTEGRAND  FOR  STRUVE  FUNCTION 

3  CKERNL  -  DSIN(PARAM*DC0S(Q))*(DSIN(Q))**2. 

RETURN 

END 

C  INTEGRATION  SUBROUTINE 

SUBROUTINE  INTGRl (XL , XU , N, NPARM , PAS S , ANS) 

REAL*8  ANS, H, XL. XU, CKERNL, GAU SS ,DFLOAT, PASS 
EXTERNAL  CKERNL, GAUSS 
ANS  -  0. 

H  -  (XU-XL)/DFLOAT(N) 

IF(H.LE.0.0  DO)  GO  TO  2 
DO  1  I-1,N 

1  ANS-ANS+GAUSS(XL+(DFL0AT(I)-1 . )*H, XU- 
1 (DFLOAT(N)-DFLOAT(I))*H,NPARM,PASS) 

RETURN 

2  ANS  -  0.0  DO 
RETURN 

END 

C  GAUSSIAN  ROUTINE 

DOUBLE  PRECISION  FUNCTION  GAUSS ( XI , X2 , MPARM, PASS ) 

REAL*8  XI ,X2,T,E,C, CKERNL, PASS, CONST 
EXTERNAL  CKERNL 

DATA  CONST/0.1127016653792583D+00/ 

T  -  X2  +  XI 
E  -  X2  -  XI 
C  -  CONST  *  E 

GAUSS-  E/18. *(5.*CKERNL(C+X1, MPARM, PASS)+8.*CKERNL(T/2. , 
1 MPARM, P AS S )+5 .*CKERNL(X2-C, MPARM, PASS) ) 

RETURN 

END 
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